TIT-HEP-519 
UCLA/04/TEP/17 
|hep-th/0405101| 



Supersymmetric Gauge Theories with Flavors 

and Matrix Models 

Changhyun Ahn 1,2 , Bo Feng 1 , Yutaka Ookouchi 3 and Masaki Shigemori 4 

School of Natural Sciences, Institute for Advanced Study, Olden Lane, Princeton NJ 08540, USA 
2 Department of Physics, Kyungpook National University, Taegu 702-701, Korea 
3 Department of Physics, Tokyo Institute of Technology, Tokyo 152-8511, Japan 
^Department of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1547, USA 



We present two results concerning the relation between poles and cuts by using the example 
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Abstract 



1 Introduction 



String theory can be a powerful tool to understand four dimensional supersymmetric gauge 
theory which exhibits rich dynamics and allows an exact analysis. In , using the generalized 
Konishi anomaly and matrix model [2], J\f — 1 supersymmetric U(N C ) gauge theory with 
matter fields in the adjoint, fundamental and anti-fundamental representations was studied. 
The resolvents in the quantum theory live on the two-sheeted Riemann surface defined by 
the matrix model curve. Their quantum behavior is characterized by the structure around 
the branch cuts and poles, which are related to the RR flux contributions in the Calabi-Yau 
geometry and flavor fields, respectively. A pole associated with flavor on the first sheet is 
related to the Higgs vacua (corresponding to classical nonzero vacuum expectation value of 
the fundamental) while a pole on the second sheet is related to the pseudo-confining vacua 
where the classically vanishing vacuum expectation value of the fundamental gets nonzero 
values due to quantum correction. 

It is known 1J that Higgs vacua and pseudo-confining vacua, which are distinct in the 
classical theory, are smoothly transformed into each other in the quantum theory. This 
transition is realized on the Riemann surface by moving poles located on the second sheet to 
pass the branch cuts and enter the first sheet. This process was analyzed in pQ at the off-shell 
level by fixing the value of glueball fields during the whole process. However, in an on-shell 
process, the position of poles and the width and position of branch cuts are correlated (when 
the flavor poles are moved, the glueball field is also changed). It was conjectured in £Q that for 
a given branch cut, there is an upper bound for the number of poles (the number of flavors) 
which can pass through the cut from the second sheet to the first sheet. 

Our first aim of this paper is to confirm this conjecture and give the corresponding upper 
bound for various gauge groups (in particular, we will concentrate on the U(N C ) gauge group). 
The main result is that if Nf > N c , the poles will not be able to pass through the cut to the 
first sheet where N c is the effective fluxes associated with the cut (and can be generalized to 
other gauge groups). 

Another important development was made in jHj, which was inspired by j3]. In which 
we will refer to as IKRSV, it was shown that, to correctly compute the prediction of string 
theory (matrix model), it is crucial to determine whether the glueball is really a good variable 
or not. A prescription was given, regarding when a glueball field corresponding to a given 
branch cut should be set to zero before extremizing the off-shell glueball superpotential. The 
discussion of IKRSV was restricted to Af = 1 gauge theories with an adjoint and no flavors, 
so the generalization to the case with fundamental flavors is obviously the next task. 
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Our second aim of this paper is to carry out this task. The main result is the following. 
Assuming Nf poles around a cut associated with gauge group U(N Cti ), when Nf > N c>i there 
are situations in which we should set Si = in matrix model computations. More concretely, 
situations with Si — belong to either of the following two branches: the baryonic branch 
for N Cj i < Nf < 2N Cj i, or the r = N c> i non-baryonic branch for Nf > 2N C ^. Moreover, when 
Si — 0, the gauge group is completely broken and there should exist some extra, charged 
massless field which is not incorporated in matrix model. 

In section EJ as background, we review basic materials for M = 1 supersymmetric U(N C ) 
gauge theory with an adjoint chiral superfield, and Nf flavors of quarks and anti-quarks. 
The chiral operators and the exact effective glueball superpotential are given. We study the 
vacuum structure of the gauge theory at classical and quantum levels. We review also the 
main results of IKRSV. In addition to all these reviews, we present our main motivations of 
this paper. 

In section |3J we apply the formula for the off-shell superpotential obtained in pQ to the 
case with quadratic tree level superpotential, and solve the equation of motion derived from 
it. We consider what happens if one moves Nf poles associated with flavors on the second 
sheet through the cut onto the first sheet, on-shell. Also, in subsection 13.41 we briefly touch 
the matter of generalizing IKRSV in the one cut model. 

In section EJ we consider cubic tree level superpotential. On the gauge theory side, the 
factorization of the Seiberg-Witten curve provides an exact superpotential. We reproduce 
this superpotential by matrix model, by extremizing the effective glueball superpotential with 
respect to glueball fields after setting the glueball field to zero when necessary. We present 
explicit results for U(3) theory with all possible breaking patterns and different number of 
flavors (Nf = 1,2,3,4, and 5). 

In section El after giving concluding remarks, we repeat the procedure we did in previous 
sections for SO(N c )/USp(2N c ) theories, briefly. 

In the appendix, we present some proofs and detailed calculations which are necessary for 
the analysis in sectional 

Since string theory results in the dual Calabi-Yau geometry are equivalent to the matrix 
model results, we refer to them synonymously through the paper. There exist many related 
works to the present paper. For a list of references, we refer the reader to [3] . 

2 Background 

In this section, we will summarize the relevant background needed for the study of M = 1 
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supersymmetric gauge theory with matter fields. 



2.1 The general picture of matrix model with flavors 

The generalized Konishi anomaly interpretation to the matrix model approach for TV = 1 
supersymmetric gauge theory with flavors was given in [HI QJ. Here we make only a brief 
summary on some points we will need. 

Let us consider Af = 1 supersymmetric U (N c ) gauge theory, coupled to an adjoint chiral 
superfield $, Nf fundamentals Q* , and Nf anti-fundamentals Qj. The tree level superpoten- 
tial is taken to be 

W tree = TrWW+^QfnffWQ* > (2- 1 ) 

/•/ 

where the function W(z) and the matrix m Jz) are polynomials 

n fe+l _ l+l 

k=0 k=l 

Classically we can have the "pseudo-confining vacua" where the vacuum expectation values of 
Q, Q are zero, or the "Higgs vacua" where the vacuum expectation values of Q, Q are nonzero 
so that the total rank of the remaining gauge groups is reduced. These two vacua, which seem 
to have a big difference classically, are not fundamentally distinguishable from each other in 
the quantum theory and in fact can be continuously transformed into each other, as we will 
review shortly, in the presence of flavors |T]. 

Supersymmetric vacua of gauge theory are characterized by the vacuum expectation values 
of chiral operators [7j . They are nicely packaged into the following functions called resolvents 
0Ej: 1 



p , s i / T w a w a \ ro _. 

1 

where W a is (the lowest component of) the field strength superfield. Classically, R(z) vanishes 
while T(z), M(z) have simple poles on the complex 2-plane at infinity and at the eigenvalues 

1 We set w a (z) = to zero because in supersymmetric vacua w a (z) = 0. 



M V) f f = (Qt—zQ f ) ( 2 - 4 ) 
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of $. Each eigenvalue of $ is equal to one of zeros of W'(z) or B(z), where 

n L 

W'{z)=g n \[{z-a i ), B(z)=detm(z)=B L l[(z-z I ). (2.5) 

i=i 1=1 

In the pseudo-confining vacuum, every eigenvalue of $ is equal to for some i. On the other 
hand, in the Higgs vacuum, some eigenvalues of $ are equal to zj for some /. 

In the quantum theory, the resolvents (|2.2j) . (J2.3|) and (|2.4j) are determined by the gener- 
alized Konishi anomaly equations [HI El U| : 

\W'(z)T{z)\- + Tr[m'(z)M(z)}_ = 2R(z)T(z), 
[W'(z)R(z)}- = R{z)\ 
i(M(z)m(z)y f ' f ]- = R(z)S f f ', (2 ' 6) 
[(m(z)M(z)f 7 U = R{z)5 J r 

where the notation [ ]_ means to drop the nonnegative powers in a Laurent expansion in z. 
From the second equation of (j2.6j) . one obtains 



R{Z) = ~ (W'(Z) - y/W'(z)* + /(*)) 



where f(z) is a polynomial of degree (n — 1) in z. This implies that in the quantum theory 
the zeros z — a$ (i — 1, 2, . . . , n) of W'(z) are blown up into cuts along intervals 2 [a~, a^] 
by the quantum effect represented by f(z), and the resolvents ()2.2j) - (j2.4|) are defined on a 
double cover of the complex z-plane branched at the roots af of W'(z) 2 + f(z). This double 
cover of the z-plane can be thought of as a Riemann surface £ described by the matrix model 
curve 

E: yl = W\zf + f{z). (2.7) 

This curve is closely related to the factorization form of M = 2 curve in the strong coupling 
analysis. 

Every point z on the z-plane is lifted to two points on the Riemann surface E which we 
denote by q and q respectively. For example, zj is lifted to qj on the first sheet and g} on the 
second sheet. We write the projection from E to the z-plane as zj = z(qi) = z(qi), following 
the notation of pQ. 

The classical singularities of the resolvents T(z), M(z) are modified in the quantum theory 

to the singularities on E, as follows. For T(z), the classical poles at zi are lifted to poles at qi 

2 af are generally complex and in such cases we take Ai to be a straight line connecting a~ and af. Note 
that there is no physical meaning to the choice of the cut; it can be any path connecting and af . 
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or qj, depending on which vacuum the theory is in, while the classical poles at <Zj with residue 
N c>i are replaced by cuts with periods § A T(z)dz = N Cji . For M(z), the classical poles at 
zj are also lifted to poles at qj or qj. More specifically, by solving the last two equations of 
(|2.6j) . one can show [T] 



^ 1 y—, (1 — rj)R(qj) 1 / dx y-^ rjR(qj) 1 / dx 

m ( Z ) jT[ ( Z ~ Z l) 27li hi m ( X ) J^_( Z ~ Zl ) 27ri hi m ( X "> ' 



where (qi,qi) are the lift of zj to the first sheet and to the second sheet of E, and n = for 
poles on the second sheet and rj = 1 for poles on the first sheet. Furthermore, for T(z), by 
solving the first equation of ()2.6|) . 

T(z) _ B'(z) A (1 - 2r J )y( gJ ) c(z) 

J ^ 2S(z) Z. 2 ^ )(z _^) + y(z) « ^ 

where 

c M = /lV ly 'W-fW\-lt ly ' W "' y ' fe) - P-10) 

\ / 2 f^ 

Practically it is hard to use ()2.1()j) to obtain c(z) and we use the following condition instead: 

^- I T{z)dz = N Cii . (2.11) 
Finally, the exact, effective glueball superpotential is given by [T] 

L /-Ao 1 L /-Ao 



'•/; 

n-1 



1 /• 1 /" 1 /" 

W cS = --J2N c J y(z)dz--"£(l-ri) y(z)dz - -"£ri 

2 i= i j bi 1 1=1 hi 1 1=1 hi 

L 

+ -{2N C - L)W(A ) + W ( Zl ) ~ ni ( 2N c ~ L)S + 2nir S + 2m 6^, (2.12) 
2 2 

7=1 8=1 

where 

2vr.ro = log j . (2.13) 

Here, A is the cut-off of the contour integrals, A is the dynamical scale, S = Y^i=i an d 
ft, G Z. B[ is the regularized contour from A to A through the z-th cut and A and A 
are the points on the first sheet and on the second sheet, respectively. The glueball field is 
defined as 



S t = ^-<( R(z)dz. 
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In the above general solutions ([2.8)1 . ()2.9|) . we have 77 = 1 or rj = 0, depending on whether 
the pole is on the first sheet or on the second sheet of £, respectively. The relation between 
these choices of 77 and the phase of the system is as follows. Let us start with all 77 = 0, i.e., 
all the poles are on the second sheet. This choice corresponds to the pseudo-confining vacua 
where the gauge group is broken as U(N C ) — > Yli=i U(N Ct i) with Yli=i N c> i = N c . Now let 
us move a single pole through, for example, the n-th cut to the first sheet. This will break 
the gauge group as niLi U(N C)i ) — ► n^i 1 U(N c< i) x U(N c ^ n — 1). Note that the rank of the 
last factor is now (iV Cjn — 1) so that ^™ =1 iV C) j = N c — 1 < iV c . Namely, the gauge group is 
Higgsed down. In this way, by passing poles through cuts, one can go continuously from the 
pseudo-confining phase to the Higgs phase, as advocated before. 

However, if we consider this process of passing poles through a cut to the first sheet on- 
shelly then there should be an obstacle at a certain point. For example, if initially we have 
N c>n = 1, after passing a pole we would end up with an U (0). This sudden jump of the number 
of £7(1) 's in the low energy gauge theory is not a smooth physical process, because the number 
of massless particles (photons) changes discontinuously. So we expect some modifications to 
the above picture. In jT], it was suggested that in an on-shell process, the n-th cut will close 
up in such a situation so that the pole cannot pass through. It is one of our motivations to 
show that this is indeed true. More precisely, the cut does not close up completely and the 
pole can go through a little bit further to the first sheet and then will be bounced back to the 
second sheet. 

2.2 The vacuum structure 

In the last subsection we saw that different distributions of poles over the first and the 
second sheets correspond to different phases of the theory. In this subsection we will try to 
understand this vacuum structure of the gauge theory at both classical and quantum levels 
for a specific model (for more details, see ^OJEIIEIEIIEIEI])- For simplicity we will focus 
on U (N c ) theory with Nf flavors and the following tree level superpotential 3 

Nf 

W trcc = ^m A Tr $ 2 - Q/($ + m^Q 1 . (2.14) 

1=1 

This corresponds to taking polynomials in (|2.1j) as 

W{z) = ^V, m 7 j(z) = -(z + m f )8j. 

3 We used the convention of JH| f° r the normalization of the second term. Different choices are related to 
each other by redefinition of Q and Q. 
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All Nf flavors have the same mass to/, and the mass function defined in ()2.5|) is given by 

B(z) = (-l) N f(z + m f ff. 
Therefore, poles associated with flavors are located at 



Zl 



-TO/ = z f , 



1,2,. ..,N f . 



(2.15) 



In the quantum theory, some of these poles are lifted to g/ on the first sheet and others are 
lifted to qf on the second sheet. 

The D- and F-flatness for the superpotential ()2.14|) is given by 



o = [$,$t], o = Qg t -g t g, o = m A $-gg, o = ($ + m,)Q = Q($ + 



Tfl 4 



Solutions are a little different for m/ ^ and rrif = 0, because m/ = is the root of 
W'(z) = z. The case of W'{— TO/) = was discussed in PUJ E2] which we will refer to as the 
classically massless case. 

In the rrif ^ case, the solution is given by 



Q 



-m f I KxK 

0(AT c _^) x (Ar c _ii-) 



Arxk 






{N c -K)x(N f -K) 



'A 



KxK 





B 







Kx(N f -K) 



(N c ~K)x(N f -K) 



(2.16) 



with 



-m f m A I KxK = A% AA^ = A^A + B* l B. 
The gauge group is Higgsed down to U(N C — K) where 

K mf¥ Q < min (JV C , Nf) . 

To understand the range of K mfJ L , first note that the $ breaks the gauge group as U(N C ) — > 
U(K) x U(N C — f^). Now the U(K) factor has effectively iV/ massless flavors and because 
(QQ) ^ 0, f/(K) is further Higgsed down to 17(0). 

For rrif = 0, we have $ = and Q,Q are still of the above form ()2.16|) with one special 
requirement: A = 0. Because of this we have 



K 



rrif=0 



< min ( N, 



N f 



where [ ] means the integer part. The integer K m =0 precisely corresponds to the r-th branch 
discussed in jl()l IT3] . The rrif — case is different from the rrif ^ case as follows. First, 
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$ does not break the U(N C ) gauge group, i.e., U(N C ) — > U(N C ). Secondly, the r-th branch is 
the intersection of the Coulomb branch in which (QQ) = and the Higgs branch in which 
(QQ) 0) whereas for rrif ^ the vacuum expectation value (QQ) must be nonzero and 
the gauge group must be Higgsed down. For these reasons, K m ,^o and K mf=0 have different 
ranges. 

The above classical classification of r-th branches is also valid in the quantum theory 
(including the baryonic branch). 

The quantum r-th branch can also be discussed by using the Seiberg-Witten curve. In 
the r-th branch, the curve factorizes as 

= (x + m f ) 2r [P 2 Nr _M ~ 4A 2Nc ~ Nf (x + m f ) N f- 2r ] . 

Because N c — r > (coming from P^ c _ r (x)) and Nf — 2r > (coming from the last term), 
we have r < N c and r < Nf/2, which leads to the range 

(2.17) 

The relation between this classification of the Seiberg-Witten curve and the above classifica- 
tion of r-branches, in the rrif ^ and rrif = cases, is as follows. In the rrif = case, we 
have one-to-one correspondence where the r is identified with K m . = Q. In the rrif ^ case, 
for a given r of the curve, there exist two cases: either K mf ^ = r for K mfJ L < [Nf/2], or 
K mf ^ = N f -r for K mf¥o > [N f /2\. 

2.3 The work of IKRSV 

Now we discuss another aspect of the model. In the above, we saw that there is a period 
condition (|2.11j) for T(z). So, if for the z-th cut we have § A T(z)dz = N C) i = 0, then it seems 
that, in the string theory realization of the gauge theory, there is no RR flux provided by 
D5-branes through this cut and the cut is closed. Because of this, it seems that we should set 
the corresponding glueball field Si — 0. Based on this naive expectation, Ref. jlj calculated 
the effective superpotential of USp(2N c ) theory with an antisymmetric tensor by the matrix 
model, which turned out to be different from the known results obtained by holomorphy and 
symmetry arguments (later Refs. [T71 ITS] confirmed this discrepancy). 

This puzzle intrigued several papers [HI 1201 El HH1 12H 123 1231 12H In particular, in [20], 
it was found that although N c ^ = 0, we cannot set = 0. The reason became clear by 
later studies. Whether a cut closes or not is related to the total RR flux which comes from 



r < min ( N c , 



■Nf_ 
2 
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both D5-branes and orientifolds. For USp(2N c ) theory with antisymmetric tensor, although 
the RR flux from D5-branes is zero, there exists RR flux coming from the orientifold with 
positive RR charges, thus the cut does not close. That the cut does not close can also be 
observed from the Seiberg-Witten curve ^3] where for such a cut, we have two single roots in 
the curve, instead of a double root. All these results were integrated in 3J for Af = 1 gauge 
theory with adjoint. Let us define N c = N c for U(N C ), N c /2 - 1 for SO(N c ) and 2N C + 2 for 
USp(2N c ). Then the conclusion of ^ can be stated as 

If N c > 0, we should include Si and extremize W e ff(Si) with respect to it. On the 
other hand, if N c < 0, we just set Si = instead. 

In it was argued that this prescription of setting Si — can be explained in string 
theory realization by considering an extra degree of freedom which corresponds to the D3- 
brane wrapping the blown up S 3 and becomes massless in the S — > limit [22] • Our second 
motivation of this paper is to generalize this conclusion to the case with flavors. We will 
discuss the precise condition when one should set Si = in order to get agreement with the 
gauge theory result, in the case with flavors. 

2.4 Prospects from the strong coupling analysis 

Before delving into detailed calculations, let us try to get some general pictures from the 
viewpoint of factorization of the Seiberg-Witten curve. Since we hope to generalize IKRSV, 
we are interested in the case where some S{ vanish. Because Si is related to the size of a 
cut in the matrix model curve, which is essentially the same as the Seiberg-Witten curve, we 
want some cuts to be closed in the Seiberg-Witten curve. Namely, we want a double root in 
the factorization of the curve, instead of two single roots. 

For U(N C ) theory with Nf flavors of the same mass rrif = —Zf, tree level superpotential 
(12.11) . and breaking pattern U(N C ) — > Yl™ =1 U(N Cti ), the factorization form of the curve is [T3] 



P Nc {zf - AA 2N «- N i{z - z f ff = H N „ n (z) 2 F 2n (z), 

F 2n (z) = W'(z) 2 + f n _ 1 (z), (2,18) 

where the degree 2n polynomial F2 n {z) = W'(z) 2 + / n _i(z) generically has 2n single roots. 
How can we have a double root instead of two single roots? 

For a given fixed mass, for example Zf = a 1; there are three cases where we have a 
double root, as follows, (a) There is no U(N Ct i) group factor associated with the root a±, 
namely N Ci i = 0. (b) The U(N Ct i) factor is in the baryonic branch. This can happen for 
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N c> i < Nf < 2N Cj i. (c) The U(N Ct i) factor is in the r-th non-baryonic branch with r = iV Cjl . 
This can happen only for Nf > 2N C ^. Among these three cases, (b) and (c) 13J are new for 
theories with flavors, and will be the focus of this paper. However, it is worth pointing out 
that the factorization form in the cases (b) and (c) are not the one given in (J2.18|) for a fixed 
mass, but the one given in (jA.3|) . 

Can we keep the factorization form (j2.18|) while having an extra double root? We can, but 
instead of a fixed mass we must let the mass "floating," which means the following. There 
will be multiple solutions to the factorization form (|2.18|) . and for any given solution the 2n 
single roots of F 2n (x), denoted by af, i = 1, . . . ,n, are functions of Zf. Now, we tune Zf so 
that af(zf) = a^(zf), i.e., so that two single roots combine into one double root. Since for 
different solutions this procedure will lead to different values of Zf, we call this situation the 
"floating" mass. 

Now we have two ways to obtain extra double roots: one is with a fixed mass, but to 
go to the baryonic or the r = iV C) i non-baryonic branch, while the other is to start with a 
general non-baryonic branch but using a floating mass. In fact it can be shown that these two 
methods are equivalent to each other when the double root is produced. In the calculations 
in section 0] we will use the floating mass to check our proposal. 

3 One cut model — quadratic tree level superpotential 

In this section we will study whether a cut closes up if one tries to pass too many poles 
through it. If the poles are near the cut, the precise form of the tree level superpotential 
(namely the polynomials W{z), f(z)) is inessential and we can simplify the problem to the 
quadratic tree level superpotential given by (|2.14|) . For this superpotential, we will compute 
the effective glueball superpotential using the formalism reviewed in the previous section. 
Then, by solving the equation of motion, we study the on-shell process of sending poles 
through the cut, and see whether the poles can pass or not. 

Also, on the way, we make an observation on the relation between the exact superpotential 
and the vacuum expectation value of the tree level superpotential. 

3.1 The off-shell W eS , M{z) and T(z) 

First, let us compute the effective glueball superpotential for the quadratic superpotential 
(j2.14j) . The matrix model curve ()2.7|) is related to W'(z) in this case as 

V 2 m = W'{zf + h{z) = m 2 A z 2 - M = m\{z 2 - /I), /Z = JL. (3.1) 
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Let us consider the case with K poles on the first sheet at qj = qj, for which rj — 1, and 
with (Nf — K) poles on the second sheet at qi = qj, for which 77 = (recall that (<?/, <?/) is 
the lift of Zf defined in (|2.15|) ). Using the curve (|3.1|) and various formulas summarized in 
the previous section, one can compute 



S 

n 



n f,i 



n f,i 



1 

2iri 



R{z)dz 



Am a ' 



y(z)dz = m A Al — 2S — 2S log 



S ' 



A 



y(z)dz — — I y(z)dz 



A 



2 

Ao 



2Slog-^ + 25 
Ao 



-i»Mi + ii/' 



4S 



y{z)dz 
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An 



y{z)dz 



-25 log -^ + 25 
A 



-log 



2 2 



1- 



4S \ 



m A zj 



m A Zj 
AS 



m A zj 
AS 



AS 

m A z'f 



+ 7: 



AS 



m A Zj 



where we dropped C(l/A ) terms. We have traded qi, qj for zj in the square roots, so that 

zj for \zi\ very 



1 45 

1 IT ~ 



the sign convention is such that 
large. Substituting this into (j2.12j) . we obtain 



25 

m A z'j 



and 



W eS (S) = S 



N c + log 



m. 



S N C 



E 5 

i>/=0 
7,77=1 



45 \ 



m A z] 



4S \ 



+ 



+ 



45 



m A z] 
AS 



AS 

m A 



AS 

m A zj 



1 

+ 2 



"4 1- 



4g 
m A z 2 j 



- 1 



(3.2) 



where A is the dynamical scale of the corresponding M = 2 gauge theory defined in (j2.13|) . 

Let us compute resolvents also. The resolvent M(z) is an Nf x Nf matrix. Using ()2.8|) . 
we find that for the J-th eigenvalue 
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Expanding Mj(z) around z = oo we can read off the following vacuum expectation values 



(QQ)i=^Y-[zi + (?ri-l) v 



(Q$Q)j = ^ [2z] + 2{2n - l)z IX /zj - f, - f, 



J2-(Q®Q-ziQQ) 



Nf7Jl A fI 



N f S. 



There is something worth noting here. One might naively expect that the exact superpotential 
is simply the vacuum expectation value of the tree level superpotential (|2.14|) . as is the case 
without flavors. However, this naive expectation is wrong! Although we have (Wtree,fund) = 
(- T,i(Q®Q ~ ziQQ)) ^ if S ^ 0, we still have 



W, 



cff,on-shell 



/ m A 
\ 2 



Tr$ z 



(3.3) 



as we will see shortly. The reason for ()3.3|) can be explained by symmetry arguments 
EZj. Although the tree level part (Wt ree ,fund) for fundamentals is generically nonzero and 
also contributes to Wiow, the contribution is precisely canceled by the dynamically generated 
superpotential W& ya , leaving only the $ part of W tree 4 - 



Let us calculate the resolvent T(z) also. For the present example, c(z) = m A (N c — 
from (|2.1()jl . Therefore, using (|2.9jl . we obtain the expansion of the resolvent T(z): 



T(z) 



— + > 
z ^ 
i 



zj 2r r - 1 

— H 

2 2 



+ 



jl(2N c - N f 



2r, - 1 



Zi<UZf-fi 



+ ■ 



(3.4) 



From this we can read off (Tr $ n ). For example, for K = we have 

<Tr$ 2 > = ^ / (z f - 



(Tt*) = %(zf-y/z)-ti J 
3.2 The on-shell solution 



H + 



(2N C - N f )Jl 



Now we can use the above off-shell expressions to find the on-shell solution. First we 
rewrite the superpotential as 



W, 



eff 



s 



N c + log 
—(Nf - K)S 



Ai 



3N C -N f ' 



-KS 



-log 



-log 

Zj__\_ 

2 2 



2 



AS 

m A 



+ 



m A z s 
AS 



z 2 - — ~- 

f m A * 



AS 
m A 



+ 



m A z f 
AS 



9 AS 

zi z f 

1 m A 



We would like to thank K. Intriligator for explaining this point to us. 



(3.5) 
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where c 1 = m^ c A 27Vc N f. We have set all z; to be 2/, i.e., all masses are the same, as 
in (j2.15|) . From this we obtain 5 



= log 



a; 



W c -Nf 



S N C 



+ Klog 



l ~2 4S 
Z f ~\ Z f ~ 



m A 



+ (N f - K) log 



r 2 _ 45 

Z f 



m A 



(3.6) 



or 



= log ( S~ Nc ) + K\og 



z f 



+ (N f - K) log 



3J-4S 



(3.7) 



where we have defined dimensionless quantities S 



and Zf 



m A A 2 4A 2 ~J A 

these massless quantities the cut is from along the interval [—2y/s, 2\/~S] 
Using (j3.ti|) or ()3.7|) it is easy to show that 



-. Note that using 



Weff. 



on-shell 



N c - 



Nf_ 
2 



m A z 2 f 



N f + (2K - N f 



1- 



AS 

rriAZ 2 f 



Also by expanding T(z) in the present case, just as we did in (j3.4jl . we can read off 

Nf + (2K - 



which gives us the relation 



m A z 2 f 




AS 

m A zj 



\ 2 



eff, on— shell 



as we promised in ()3.3|) . 

Equation (|3.6|) is hard to solve. But if we want just to discuss whether the cut closes up 
when we bring Zf — > 0, we can set K = 7 , for which (|3.6J) reduces to 



z f 



Nf 



S N f +u^ r f S 



• (3-8) 

5 It is easy to check that the equation Ij3.6|l with parameters (N c , Nf, K) is the same as the one with 
parameters (N c — r, Nf — 2r, K — r). Also from the expression H^.4(l it is straightforward to see we have 
(Tr<i> 2 ) JV Nf K = (Tr^ 2 )^ _ r N _ 2r K _ r + TZj. All of these facts are the result of the "addition map" 
observed in Furthermore, one can show that both i|3.6|) and (|3.4|l for K — are exactly the same as the 
one given by the strong coupling analysis in and the weak coupling analysis in |281 129| . 

6 As we mentioned before, for to/ ^ the allowed Higgs branch requires K < Nf but the strong coupling 
analysis gives K < Nf/2. The resolution for that puzzle is that if K > Nt/2, it is given by (Nf — K)-th 
branch of the curve. Since the same r-th branch of the curve gives both r and (Nf — r) Higgs branches, we 
expect that r and (Nf — r) Higgs branches are related. This relation is given by S = Sb 2 , Zf = Zfb with 
b = T "g A . It can be shown that with the above relation, the equation of motion of S for the K-th branch is 
changed to the equation of motion of S for the (Nf — K)-th branch. 

7 If K ^ 0, we will have U(N C ) -> U(K) x U(N C - K) and the problem reduces to of U(N C - K) with 
(Nf - 2K) flavors in the 0-th branch. 
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Here, u Nf is the Nf-th root of unity, cu Nf = e 2m / N f ^ and r = 0, 1, ... , (Nf — 1) corresponds 
to different branches of solutions. It is also amusing to note that above solution has the 
Seiberg duality [3U] where electric theory with (N c , Nf) is mapped to a magnetic theory with 
(Nf — N c ,Nf). 

With these preparations, we can start to discuss the on-shell process of passing poles 
through the cut from the second sheet. 

3.3 Passing poles through a branch cut 

Consider moving Nf poles on top of each other at infinity on the second sheet toward 
the cut along a line passing through the origin and making an angle of 9 with the real axis. 
Namely, take 8 

z(q f ) = z(q f )=pe w , p,6eR, (3.9) 

and change p from p = oo to p = — oo (see Fig.Q). This equation (J3.9)) needs some explanation. 
Remember that Zf = z(qf) = z(qf) denotes the projection from the Riemann surface £ (Eq. 
(|3.1|) ) to the z-plane. For each point z on the z-plane, there are two corresponding points: q 
on the first sheet and q on the second sheet. Although we are starting with poles at qf on the 
second sheet, we do not know in advance if the poles will pass through the cut and end up 
on the first sheet, or it will remain on the second sheet. Therefore we cannot specify which 
sheet the poles are on, and that is why we used z(qf), z(q/) in (|3.9jh instead of qf or qf. 




Figure 1: A process in which Nf poles at <jf on the second sheet far away from the cut 
approach the branch cut on the double sheeted z-plane, along a line which goes through the 
origin and makes an angle 9 with the real z axis. The "x " with dotted lines denotes the poles 
on the second sheet, moving in the direction of the arrow. The two branch points ±a are 
connected by the branch cut, which is denoted by a zigzag. 

8 Throughout this subsection, we will use the dimensionless quantities Zf, S, etc. and omit the hats on 
them to avoid clutter, unless otherwise mentioned. 
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Below, we study the solution to the equation of motion ([3.8)1 . changing p e M. from p = oo 
to p = -oo. By redefining z f , S by ^ -> Zfe ^r/(N f -2N c ) ^ 5 -> Se^Wf-™*), we can bring 
(|3.8|) to the following form: 



Zf = pe w = S* + S 1 - 1 , (3.10) 



where 



N c 

-t. 



Henceforth we will use (13.10)) . Because Zf as well as 5 is complex, the position of the branch 
points (namely, the ends of the cut), ±a, where 



a = V4S 

is also complex, which means that in general the cut makes some finite angle with the real 
shown in Fig. 

• N f = N c 

As the simplest example, let us first consider the Nf = N c (i.e., t — 1) case. We will see 
that the poles barely pass through the cut but get soon bounced back to the second sheet. 
The equation of motion ()3.1()j) is, in this case, 

Zf =pe ie = S + 1. (3.11) 

Therefore, as we change p, the position of the branch points changes according to 



AS = 2^z } - 1 = 2y/pe ie - 1. (3.12) 

Let us look closely at the process, step by step. The point is that transition between the 
first and the second sheet can happen only when the cut becomes parallel to the incident 
direction of the poles, or when the poles pass through the origin. 

(1) p ~ +00, on the second sheet: 

In this case, we can approximate the right hand side of (|3.12j) as 



I _ I M _i„-i„-»fl _ I _ip-i cose i(| + ip-i sine ) 



a = 2p*e^(l-p- 1 e- w )* ~2p*e*e-* p e =2p^e~^ p cosy e t( ^ 

Therefore, when the poles are far away, the angle between the cut and the real axis is 
approximately | > (we assume < 9 < |). Furthermore, as the poles approach (p 
becomes smaller), the cut shrinks (because of pa) and rotates counterclockwise (because 
of e^ p lsm6) ). This corresponds to Fig. |2K 
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(2) Because the cut is rotating counterclockwise, as the poles approach, the cut will even- 
tually become parallel to the incident direction, at some point. This happens when 

a 2 = 4{pe i9 - 1) = 4 [(p cos Q - 1) + ip sin 9} oc e 2i9 = cos 29 + i sin 29. 

By simple algebra, one obtains 

p = 2cos9, a = 2e ie . (3.13) 

Note that this is the only solution; the cut becomes parallel to the incident direction 
only once. Because < p < \a\ = 2 (we are assuming < 9 < 7r/2), by the time the 
cut becomes parallel to the incident direction, the poles have come inside of the interval 
[— 2e l9 , 2e l9 }, along which the cut extends when it is parallel to the incident direction. 
This implies that the poles cross 9 the cut at this point, and enter into the first sheet. 
Fig. |2b shows the situation when this transition is about to happen. In Fig. |2fc, the poles 
are just crossing the cut. Fig. EJi corresponds to the situation just after the transition 
happened; the poles have passed the cut and are now proceeding on the first sheet. 

Here we implicitly assumed that the cut is still rotating counterclockwise with a finite 
angular velocity, but this can be shown by expanding a around ()3.13|) as p = 2 cos 9+Ap. 
A short computation shows 

a ~ 2e^ Ap cos9 e i9 -^ A P sind 

which implies that the cut shrinks and rotates counterclockwise if we move the poles to 
the left (Ap decreases). 

(3) p ~ 0: 

If the poles proceed on the real line further, it eventually reaches the origin p = 0. By 
expanding (|3.12|) around p = 0, one obtains 

a = 2(e™ + pe i9 )^ ~ 2e -^ cos V ( 5"^ sine) . (3.14) 

Therefore the cut has a finite size (\a\ = 2) at p = and along the imaginary axis, still 
rotating counterclockwise, but now expanding. Because the cut goes through the origin, 
the poles pass through the cut again and comes back onto the second sheet (Fig. |2t)- 

(4) p ~ — oo, on the second sheet: 

If the poles have gone far past the cut so that p < 0, \p\ ^> 1, we can approximate 
9 As mentioned in footnote [3 there is no real physical meaning to the position of the cut itself. 
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Figure 2: Six configurations of the branch cut and the poles. The poles are depicted by "x " 
and moving along a line at an angle 9 with the real axis, as the arrow on it indicates. The "x " 
in solid (dotted) lines denotes poles on the first (second) sheet. The branch cut is rotating 
counterclockwise (as the arrows on its sides indicate), changing its length. 

()3.12|) . as before, as 

a = 2\p\h^\l - V -\- %e )^ ~ 2H^l p |- w ^(i+HW- 1 ^). 

Therefore, as the poles go away, the cut expands and rotates counterclockwise. The 
angle between the cut and the real axis asymptotes to (| + |) (Fig. |2f). 

In the above we assumed that < 9 < |. If | < 9 < ir, the only difference is that 
the order of steps (2) and (3) are exchanged. If 9 < 0, the cut rotates clockwise instead of 
counterclockwise. 

When is the cut shortest in this whole process? From (J3.12)) . one easily obtains 

\ a \ = 2[{p-cos9) 2 + sin 2 #] 1 / 4 > 2|sinfl| 1 / 2 . (3.15) 

Therefore, when the poles are at zj = pe td = cos9e t0 , which is between the steps (2) and 
(3) above, the cut becomes shortest. In particular, in the limit 9 — >• or 9 — > ±tt, the 
cut completely closes up instantaneously. These correspond to configurations with either a 
horizontal cut with poles colliding sideways, or a vertical cut with poles colliding from right 
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above or from right below. Actually the existence of the 5 = solution is easy to see in 
(prrrjl : it is just Zf = l,S = 0. 

Summary: for Nf = N c , when one moves poles on the second sheet from infinity along a 
line toward a cut, poles pass through the cut onto the first sheet and move away from the cut 
by a short distance. Then poles are bounced back to the second sheet again. Therefore, one 
can never move poles far away from the cut on the first sheet. During the process, in certain 
situations, the cut completely closes up. 

• N f ?N c 

Now let us consider a more general case with Nf 7^ N c . We again consider a situation 
where poles on the second sheet approach a cut. This time we will be brief and sketchy, 
because a detailed analysis such as the one we did for the Nf = N c case would be rather 
lengthy due to the existence of multiple branches, and would not be very illuminating. 

First, let us ask how we can see whether poles are on the first sheet or on the second sheet, 
from the behavior of S versus p. Because this is not apparent in the equation of motion of 
the form ()3.10|) . let us go back to 



o^^ + ^f^V^ =* s'^^pl. „ 16) 

which led to the equation (j3.1Uj) . Here the "— " ("+") sign corresponds to qf on the first 
(second) sheet. For \zf\ 2 ^> \4:S\, the square root can be approximated as J z 2 — 4S = 
Zf(l — AS/ z 2 ) 1 / 2 ~ Zf(l — 2S/z 2 ) (our sign convention was discussed above ([3.2)1 ). Therefore 
(J3.16)) is, on the first sheet, 

, z f - z f (l - 2S/z 2 f ) S , , 

2 Zf 

while on the second sheet 

Q t Zf + Zf(l-2S/z 2 ) 

S ~ — ~ Zf =^ Zf ~ Sr. (3.18) 

Now let us solve ()3.1()j) for \z/\ 3> 1. By carefully comparing the magnitude of the two 
terms in (|3.10|) . one obtains 

N f <N C (Kt) -> 



N c <N f <2N c (|<*<1) 
2N C <N f (0 < t < ~) 



Zf ~ S t 


- \s\ 




\s\ » 1 


Z f ~ S 1 ^ - 


-> \s\ 




\s\ < 1 


Zf ~ S l 


-» 151 


~ \p\T, 


\s\ » 1 


Zf ~ S 1 - 1 


- \s\ 


- bl 1 ^, 


|5| > 1 



(3.19) 
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It is easy to show that \zf\ 2 ^> \4S\ in all cases. So by using ()3.17|) . ()3.18|) . we conclude 
that the first and the third lines in (j3.19|) correspond to poles on the second sheet, while the 
second and the last lines correspond to poles on the first sheet. This implies that, only for 
Nf < N c , poles on the second sheet can pass through the cut all the way and go infinitely 
far away on the first sheet from a cut, as we will see explicitly in the examples below. For 
N c < Nf < 2N C , if one tries to pass poles through a cut, then either poles will be bounced 
back to the second sheet, or the cut closes up before the poles reach it. For Nf > 2N C , there 
is no solution corresponding to poles moving toward a cut from infinity on the second sheet. 
This should be related to the fact that in this case glueball S is not a good IR field. Therefore, 
this one cut model is not applicable for Nf > 2N C . 

To argue that these statements are true, rather than doing an analysis similar to the one 
we did for the Nf = N c case, we will present some explicit solutions for some specific values 
of t = jfj and 9, and argue general features. Before looking at explicit solutions, note that 
there are multiple solutions to Eq. ()3.10j) which can be written as 

Zf = (S^) Nc + (S^) N f~ Nc . (3.20) 

From the degree of this equation, one sees that the number of the solutions to (|3.20j) is: 

Nf < N c 2N C - Nf solutions, 

N c < N f < 2N C N c solutions, 

2N C < N f N f - N c solutions. 

Therefore, if we solve the equation of motion ([3.10jh in general we expect multiple branches 
of solutions. We do not have to consider the Nf branches of the root S N f , because it is taken 
care of by the phase rotation we did above (|3.1U|) . 

Now let us look at explicit solutions for U{2) example. For Nf = \N C (t = 2), the solution 
to the equation of motion (|3.1U|) is 

27+ J 729 - IO84 \ /27 + «/729- 1082}\ 

5 J "I 5 ) • 

where three branches of the cubic root are implied. In Fig. El we plotted \S\ versus p for these 
branches, for a randomly chosen value of the angle of incidence, 9 = 7r/6. Even if one changes 
9, there are always three branches whose general shapes are similar to the ones in Fig. El 
These three branches changes into one another when 9 is changed by 27r/3. One can easily 
see which branch corresponds to what kind of processes, by the fact that on the first sheet 
|tS| C 1 as \p\ — > oo, while on the second sheet \S\ > 1 as \p\ — > oo. The three branches 
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correspond to the process in which: i) poles go from the second sheet to the first sheet through 
the cut, without any obstruction, ii) poles go from the first sheet to the second sheet (this is 
not the process we are interested in), and iii) poles coming from the second sheet get reflected 
back to the second sheet. Note that the cut has never closed in all cases, because \S\ is always 
nonvanishing. 




branch i) branch ii) branch iii) 



Figure 3: The graph of \S\ versus p for Nf = |iV c (t = 2), 9 = tt/6. The vertical axis is \S\ 
and the horizontal axis is p. Although we are showing just the 9 = tt/6 case, there are similar 
looking three branches for any value of 9, which change into one another when 9 is changed 
by 2tt/3. 



Similarly, for Nf = |iV c (t = |), the solution to the equation of motion (|3.10|) is 

S=±[-3z,-l±{Zf + l)y/tePn]. 

This time there are two branches, which change into each other when the 9 is changed by n. 
We plotted \S\ versus p for 9 = n/2 in Fig. 0] It shows two possibilities: i) poles coming from 
the second sheet get reflected back to the second sheet, for which \S\ ^ as p — > 0, ii) the 
cut closes up before poles passes through it, for which \S\ — > as p — > 0. 




-1 -0.7^0.5-0.25 0.25 0.5 0.75 1 -1 -0.75-0.5-0.25 0.25 O.o 0.75 1 



branch i) branch ii) 

Figure 4: The graph of \S\ versus p for N f = §iV c (t = \) or N f = 3N C (t = ±), for 9 = vr/2. 
The two values oft give the same graph. The vertical axis is \S\ and the horizontal axis is 
p. Although we are showing just the 9 = tt/2 case, there are similar looking two branches for 
any value of 9, which change into each other when 9 is changed by it, 
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The Nf = 3N C (t = g) case is also described by the same Fig.HJ However, as we discussed 
below (|3.19p . it does not correspond to a process of poles approaching the cut from infinity 
on the second sheet; it corresponds to poles on the first sheet and we cannot give any physical 
interpretation to it. 

These demonstrate the following general features: 

• For Nf < N c , one can move poles at infinity on the second sheet through a cut all the 
way to infinity on the first sheet without obstruction, if one chooses the incident angle 
appropriately. If the angle is not chosen appropriately, the poles will be bounced back 
to the second sheet. 

• For N c < Nf < 2N C , one cannot move poles at infinity on the second sheet through a 
cut all the way to infinity on first sheet. If one tries to, either i) the cut rotates and 
sends the poles back to the second sheet, or ii) the cut closes up before the poles reach 
it. 

• For 2N C < Nf, the one cut model does not apply directly. 

Nf = 2N C is an exceptional case, for which the equation of motion (j3.10|) becomes 

z f = 2S 1/2 . (3.21) 
Therefore \S\ ^ as Zf — > 0, and the cut always closes before the poles reach it. 
Subtlety in S = solutions 

If p — 0, or equivalently if Zf = 0, there is a subtle, but important point we overlooked in the 
above arguments. For Zf = 0, the superpotential ()3.5|) becomes 



W, 



eff 



s 



s 



N f \ , f (-l) N f/ 2 mf- Nf/2 A 2N °- N f' 

^-^j+iog r s i_ Nf/2 



N r 



N, 



log 



'/ -.\-jW2 N c -N f /2 . 2N C -N f ' 

( — 1) f> m A 1 A 1 



1 + l0Kl ^ 



S N c -N f /2 

+ 2mToS 



+ 2niT S 



(3.22) 



with A 



3(N c -N f /2) 



(—1) N f^ 2 m^ c N f' 2 f\ 2Nc N f . Only from here to (j3.23|) . S means the 
dimensionful quantity (S = m A A 2 S; see below (|3.7jl ). In addition, in the second line of (|3.22|) . 
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we rewrite the renormalized scale A in terms of the bare scale A and the bare coupling To 
using the relation (j2.13j) . In our case, = (—l) N f. The equation of motion derived from 



[no solution Nf = 2N C . 

For Nf < N c , (j3.23|) is consistent with the fact that the \S\ versus p graphs in Fig. El all 
go through the point (p, \ S\) = (0, 1) (now 5* means the dimensionless quantity). Also for 
Nf = N c , ()3.23|) is consistent with the result (|3.14J) (|a| = 2, so \S\ = 1). On the other hand, 
for Nf > N c , ()3.23|) implies that we should exclude the origin (p, \S\) = (0,0) from the \S\-p 
graphs in Fig. EJ which is the only S = solution (this includes the Nf = 2N C case (j3.21j) ). 10 

Therefore, the above analysis seems to indicate that, for Nf > N c , the S = solution at 
p = 0, or equivalently Zf = is an exceptional case and should be excluded. On the other 
hand, as can be checked easily, gauge theory analysis based on the factorization method shows 
that there is an S = solution in the baryonic branch. Thus we face the problem of whether 
the baryonic 5* = branch for Nf > N c can be described in matrix model, as alluded to in 
the previous discussions. 

Note that, there is also an S = solution for N c = Nf in certain situations, as discussed 
below (j3.15j) . For this solution, which is in the non-baryonic branch, there is no subtlety in 
the equation of motion such as ()3.23|) . and it appears to be a real on-shell solution. This will 
be discussed further below. 

3.4 Generalization of IKRSV for the one cut model 

In the above and in subsection 12.41 we argued that for Nf > N c the 5* = solutions are 
real, on-shell solutions based on the factorization analysis. More accurately, there are two 
cases with S = 0: the one in the maximal non-baryonic branch with Nf > 2N C and the other 
one in the baryonic branch with N c < Nf < 2N C . The case of non-baryonic branch cannot 
be discussed in the one cut model, which is applicable only to Nf < 2N C . On the other 

10 One may think that if one uses the first line of $7Zty . then for Nf = 2N C , W cB = 51og[(-l)^ /2 ] and 
there are solutions for some Nf. However, the glueball superpotential that string theory predicts \'62\ is 
the third line of (|3.22() which is in terms of the bare quantities Ao and To. If Nf = 2N C , then the log term 
vanishes and one cannot define a new scale A as we did in H2.13JI to absorb the linear term 2-kitqS. 




and the solution is 




2- 



(3.23) 
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hand, the baryonic one did show up in the previous subsection, but we just saw above that 
those solutions should be excluded by the matrix model analysis. What is happening? Is it 
impossible to describe the baryonic branch in matrix model? 

Recall that the glueball field S has to do with the strongly coupled dynamics of U(N C ) 
theory. That S = in those solutions means that there is no strongly coupled dynamics any 
more, namely the U(N C ) group has broken down completely. The only mechanism for that to 
happen is by condensation of a massless charged particle which makes the U(l) photon of the 
U(N C ) group massive. Therefore, in order to make S = a solution, we should incorporate 
such an extra massless degree of freedom, which is clearly missing in the description of the 
system in terms only of the glueball S. This extra degree of freedom should exist even in the 
Nj = N c case where S = really is an on-shell solution as discussed below l)3.15J) ; we just 
could not directly see the degree of freedom in this case. 

The analysis of [3 hints on what this extra massless degree of freedom should be in the 
matrix model / string theory context. Note that, the superpotential (|3.22j) is of exactly the 
same form as Eq. (4.5) of 0, if we interpret N c — Nf/2 = N as the amount of the net RR 
3-form fluxes. In [3] it was argued that, if the net RR flux iV vanishes, one should take into 
account an extra degree of freedom corresponding to D3-branes wrapping the blown up S 3 in 
the Calabi-Yau geometry j2S] , arid condensation of this extra degree of freedom indeed makes 
S = a solution to the equation of motion. The form of the superpotential 1)3.22)1 strongly 
suggests that the same mechanism is at work for Nf = 2N C in the r = N c non-baryonic 
branch; condensation of the D3-brane makes S = a solution. Furthermore, as discussed in 

for Nf > 2N C the glueball S is not a good variable and should be set to zero. A concise 
way of summarizing this conclusion is: if the generalized dual Coxeter number h = N c — Nf/2 
is zero or negative, we should set S to zero in the r = N c non-baryonic branch. 

However, this is not the whole story, as we have discussed in subsection 12.41 As we saw 
above, we need some extra physics also for N c < Nf < 2N C in order to explain the matrix 
model result in the baryonic branch. We argue below that this extra degree of freedom at 
least in the N c < Nf < 2N C case should also be the D3-brane wrapping S 3 which shrinks to 
zero when the glueball goes to zero: S — > 0. 

The original argument of [3| is not directly applicable for Nf < 2N C because there are 
nonzero RR fluxes penetrating such a D3-brane (N ^ 0). These RR fluxes induce fundamental 
string charge on the D3-brane. Because the D3-brane is compact, there is no place for the flux 
to end on (note that this flux is not the RR one but the one associated with the fundamental 
string charge). Hence it should emanate some number of fundamental strings. If there are 
no flavors, there is no place for such fundamental strings to end on, so they should extend 
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to infinity. This fact led to the conclusion of [3] that the D3-brane wrapping S 3 is infinitely 
massive and not relevant unless N = 0. 

However, in our situation, there are places for the fundamental strings to end on - 
noncompact D5-branes which give rise to flavors pHl 132] . In particular, precisely in the 
Zf = case, where we have S = solutions for N c < Nf < 2N C , the D3-brane wrapping 
S 3 intersects the noncompact D5-branes in the S — > limit, hence the 3-5 strings stretching 
between them are massless. Therefore the D3-brane with these fundamental strings on it is 
massless and should be included in the low energy description. It is well known [33] that 
such a D-brane with fundamental strings ending on it can be interpreted as baryons in gauge 
theory. 11 Condensation of this baryon degree of freedom should make S = a solution, 
making the photon massive and breaking the U(N C ) down to U(0). The precise form of 
the superpotential for this extra degree of freedom must be more complicated than the one 
proposed in ,3J for the case without flavors. 

All these analyses tell us the following prescription: 

Using the floating mass condition that all Nf poles are on top of one branch 

point 12 on the Riemann surface, we will have an S = solution for Nf > 2N C . 

For N c < Nf < 2N C there are two solutions: one with S = in the baryonic ^ ^ 

branch and one with S ^ in the non-baryonic branch. In multi-cut cases, 

this applies to each cut by replacing N c , S with the corresponding N c j, Si for 

the cut. 

In the next section we will discuss the condition we have used in above prescription. Also 
by explicit examples, we will demonstrate that when the gauge theory has a solution with 
closed cuts (Si = 0), one can reproduce its superpotential in matrix model by setting the 
corresponding glueballs Si to zero by hand. 

4 Two cut model — cubic tree level superpotential 

Now, let us move on to U(N C ) theory with cubic tree level superpotential, where we have 
two cuts. We will demonstrate that for each closed cut we can set S = by hand to reproduce 
the correct gauge theory superpotential using matrix model. 

11 That the D3-brane wrapping S 3 cannot exist for Nf < N c can probably be explained along the same line 
as by showing that those 3-5 strings are fermionic. Also, note that the gauge group here is U(N C ), not 
SU(N C ) as in hence the "baryon" is charged under the U(l). 

12 This condition will not work for the Nf = N c> i case. 
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Specifically, we take the tree level superpotential to be 



tree 



1=1 



(4.1) 



W{z) 




m 



9 



) 



g{z - a x ){z - a 2 ). 



Here we wrote down W(z) in terms of g 2 = g, g\ = m for definiteness, but mostly we will 
work with the last expression in terms of g, a 12 . The general breaking pattern in the pseudo- 
confining phase is U(N C ) — > U(N Cjl ) x U(N C>2 ), N cA + N Cj2 = N c , N C;i > 0. In the quantum 
theory, the critical points at ai and a 2 blow up into cuts along the intervals [a]", a] 1 "] and 
[a^,a^], respectively. Namely, we end up with the matrix model curve (|2.7j) . which in this 
case is 



We will call the cuts along [a7/,a+] and [a^aj] respectively the "first cut" and the "second 
cut" henceforth. One important difference from the quadratic case is that, we can study a 
process where Nf > 2N C ^ flavor poles are near the i-th cut in the cubic case. 

As we have mentioned, our concern is whether the cut is closed or not. Also from the 
experiences in the factorization it can be seen that for Nf > N c<i , when closed cut is produced, 
the closed cut and the poles are on top of each other 13 . With all these considerations we take 
the following condition to constrain the position of the poles: 14 



If there are Si = solutions in which the closed cut and the poles are on top of each other, 
then all such solutions can be found by solving the factorization problem under the constraint 
(J4.3|) . since for such solutions z = a{ = a\ obviously. One could impose a further condition 
S\ = 0, or equivalently a x = af if one wants just closed cut solutions, but we would like to 
know that there also are solutions with Si ^ for Nf < 2iV Cjl , so we do not do that. 

To summarize, what we are going to do below is: first we explicitly solve the factorization 
problem under the constraint (|4.H|) . and confirm that the Si — solution exists when Nf > 
N c> i. Then, we reproduce the gauge theory superpotential in matrix model by setting Si = 
by hand. 

13 We do not discuss the Nf = N c ^ case where closed cut and poles are not at the same point. However 
because the 5 = solution in this case is an on-shell solution, we can reproduce the gauge theory result in 
matrix model without setting S = by hand. 

14 We could choose Zf = a± or Zf = instead of l|4.3|) . but the result should be all the same, so we take 
H4.3JI without loss of generality. 



Vm = W(z) 2 + fi(z) = g 2 {z - oj )0 - af){z - a 2 ){z - a£). 



(4.2) 



z f = a 1 . 



(4.3) 
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Before plunging into that, we must discuss one aspect of the constraint ()4.3|) and the 
r-branches, in order to understand the result of the factorization method. If one solves the 
factorization equation for a given flavor mass mf = —Zf (without imposing the constraint 
(|4.3j) ). then in general one will find multiple r-branches labeled by an integer K with range 
< K < min(iV c , [^-]) (see Eq. ([2.17]) ). This is related to the fact that the factorization 
method cannot distinguish between the poles on the first sheet and the ones on the second 
sheet. The r-branch labeled by K corresponds to distributing Nf — K poles on the second 
sheet and K poles on the first sheet. We are not interested in such configurations; we want to 
put Nf poles at the same point on the same sheet. However, as we discuss now, we actually 
do not have to worry about the r-branches under the constraint (|4.3|) . 

The r-branches with different K are different vacua in general. However, under the con- 
straint (|4.3jl these r-branches become all identical because at the branch point z = af there 
is no distinction between the first and second sheets. This can be easily seen in the matrix 
model approach. From the equation (j2.12j) . the effective glueball superpotential for the two 
cut model with Nf — K poles at q~f on the second sheet and K poles at qf on the first sheet is 

W eS = —(N^U, + iV Ci2 n 2 ) - ±(N f - K)uf - l -KTlf + 1 -{2N C - N f )W(A ) 
+ -N f W(q) - m{2N c - N f )S + 2mT S + 2mb l S l , 
where the periods are defined by 



1 



St = — I R(z)dz, Ik = 2 / y(z)dz, 



Ao 



nf = [ y(z)dz, nf 

Jq f 



Mo 


pgr Mo 


/ y{z)dz = 


/ + / 


'If 


Jq f Jq f 



y(z)dz = Aii/ + ny 



(2) 



with i = 1,2. The periods YlP , nf are associated with the poles on the first sheet and the 

(2) 

ones on the second sheet, respectively. The contour Ci for is totally on the second sheet, 

while the contour C\ for nf is from qf on the first sheet, through a cut, to A on the second 

sheet. These contours are shown in Fig. El This r-branch with K poles on the first sheet can 

be reached by first starting from the pseudo-confining phase with all Nf poles at qf (K = 0) 

and then moving K poles through the cut to qf. The path along which the poles are moved 

in this process is the difference in the contours, C\ — Gi = AC 15 . 

15 There is ambiguity in taking AC; for example we can take AC to go around af in Fig. However, 
the difference in J AC y(z)dz for such different choices of AC is 2ninSi, n e Z, which can be absorbed in 
redefinition of the theta angle and is immaterial. 
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K poles on Nf-K poles on 
the first sheet the second sheet 



Figure 5: Contours G\ and C 2 defining TLj and Ily , respectively. The part of a contour on 
the first sheet is drawn in a solid line, while the part on the second sheet is drawn in a dashed 
line. The Nf — K poles on the second sheet and the K poles on the first sheet are actually 
on top of each other (more precisely, their projections to the z-plane are.) 

When we impose the constraint (|4.3|) . then the difference AC vanishes (Fig.EJ). Therefore 
there is no distinction between C\, C2 and hence TVp = TL^p for any K. In other words, all 
K-th branches collapse 16 to the same branch under the constraint ([4.3)1 . 




Figure 6: Under the constraint ()4.3|) . contours C\ and C 2 become degenerate: C\ = C 2 . 



Now, let us explicitly solve the factorization problem under the constraint ()4.3jl . and check 
that the S = solutions exist as advertised before. In solving the factorization problem, we 
do not have to worry about the r-branches because there is no distinction among them under 
the constraint (|4.3j) . Then, we compute the exact superpotential using the data from the 
factorization and reproduce it in matrix model by setting S = by hand when S = on the 
gauge theory side. 

For simplicity and definiteness, we consider the case of U(3) gauge group henceforth. 
We consider Nf < 2N C flavors, namely 1 < Nf < 5, because Nf > 2N C cases are not 
asymptotically free and cannot be treated in the framework of the Seiberg-Witten theory. 

4.1 Gauge theory computation of superpotential 

In this subsection, we solve the factorization equation under the constraint (|4.3|) and 

compute the exact superpotential, for the system (j4.1|) with Z7 (3) gauge group and with 

16 In fact this collapse was observed in f° r SO(N c ) and USp(2N c ) gauge groups with massless flavors. 

We have seen that there are only two branches, i.e., Special branch and Chebyshev branch, which correspond 
to the baryonic branch and the non-baryonic branch in U(N C ) case. 
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various breaking patterns. 



4.1.1 Setup 

The factorization equation for U(3) theory with Nf flavors with mass rrif — —Zf is given 

by 17 dm 

P 3 (i) 2 - AA 6 - N f(z- z f ff = H.il) 2 [Wil) 2 + f x {z)] 

= H x {zf(z- ai)(z- 4){z- az)(z- a+), (4.4) 

where we set g — 1 for simplicity and W'(z) is given by (j4.1j) . The breaking pattern is assumed 
to be U(3) — > U(N c> i) x U(N Cj2 ) with N c ^ > 0. Here we used new notations to clarify the 
shift of the coordinate below. For quantities after the shift, we use letters without tildes. 
Enforcing the constraint ()4.3|) and shifting z as z = z + aj", we can rewrite this relation as 
follows: 

P 3 (z) 2 -AA 6 - N fz N f =H 1 (z) 2 p(z) 2 + A(z)] =H 1 (z) 2 [z(z - a+)(z - aj)(* - o+)] 

= H x {zf [z{z 3 + Bz 2 + Cz + D)] . (4.5) 

Because of the shift, the polynomials Ps(z) and H\(z) are different in form from Pz(z) and 
Hi(z) in (|4.4jl . We parametrize the polynomials P%{z) and H\{z) as 

P 3 (z) = z 3 + az 2 + bz + c, #i(z) = z-A (4.6) 

The parameters B, C can be written in terms of the parameters in 

W'(z) = W'(z) = z(z + m A ) = (z + a x )(z + a{ + m A ) = (z - a x ){z - a 2 ) 

by comparing the coefficients in ()4.2|) : 

-BtV3B 2 -8C -B± V3B 2 -8C 
ai- , a 2 - , 

A 2 = (2a - a x ) 2 = (a 2 - a,) 2 = 3B2 ~ 8C , (4 . 7 ) 

The ambiguity in signs in front of the square roots can be fixed by assuming a\ < a 2 . Finally 
we undo the shift by noting that 

W(z) = \l z + ^-i 2 = Z — - {a x + a 2 ) Z — + a x a 2 z + ha* - 3aja 2 ). (4.8) 
6 16 2 o 



17 From the result of Appendix i|A.3|) . the matrix model curve with flavors does not change even for Nf > N c , 
contrary to [7] . 
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N f 


breaking pattern 

T T / A T \ T T / A T \ 

U(N c j) x £7(iV c , 2 ) 


N f < iV C)1 


iV C)1 < N f < 2N cA 


2iV C)1 < Nf 


the first 
cut is 


1 


£7(2) x £7(1) 





- 


- 


open 


£7(1) x £7(2) 





- 


- 


open 


2 


£7(1) x £7(2) 


- 


- 





closed 


£7(2) x £7(1) 





- 


- 


open 


3 


£7(1) x £7(2) 


- 


- 





closed 


£7(2) x £7(1) 


- 





- 


closed 


£7(2) x £7(1) 









open 


4 


£7(1) x £7(2) 






o 


closed 


£7(2) x £7(1) 









closed 


5 


£7(1) x £7(2) 









closed 


£7(2) x £7(1) 









closed 



Table 1: The result of factorization of curves for £7(3) with up to Nf = 5 flavors. "O 
denotes which inequality Nf and N Cj ± satisfy. 



With all this setup, we can compute the superpotential as follows. First we factorize the 
curve according to (|4.5jl . Then we find the Casimirs £7i,£7 2 ,£73 from P^{z) 18 and solve for 
ai, ci2 using the last equation of (|4.7|) . Finally we put all these quantities into ()4.8j) to get the 
effective action as 

W low = (Tr W($)) = £7 3 - (ax + a 2 )£7 2 + a 1 a 2 £7 1 + -^(a\ - Za\a 2 ). (4.9) 

o 

Here, for a±, a 2 , one can use the first two equations of ()4.7jl . We also want to know whether 
the first cut (the one along the interval [o7/, a^}) is closed or not; we expect that the cut closes 
if we try to bring too many poles near the cut. As is obvious from (J4.5)) . this can be seen 
from the value of D. If D = 0, the cut is closed, while if D ^ 0, the cut is open. 

4.1.2 The result of factorization problem 

We explicitly solved the factorization problem for £7(3) gauge theories with Nf = 1, 2, . . . , 5 
and summarized the result in Table 1. Let us explain about the table. £7(iV Cjl ) x U(N C ^) 
denotes the breaking pattern of the £7(3) gauge group. The hat on the first factor means that 

18 From the coefficients of Pz{z), namely a, b and c, one can compute the Casimirs Uk = r (Tr$ fe ) using 
the quantum modified Newton relation Appendix (|A.lf) as explained in Appendix IO 



29 



the pole is at one of the branch points of the first cut (Eq. ([4.3)1 ) which is associated with 
the first factor U(N C) x). Of course this choice is arbitrary and we may as well choose U(N Cj2 ), 
ending up with the same result. Finally, whether the cut is closed or not depends on whether 
D = or not, as explained below (|4.9|) . 

Now let us look carefully at Table 1, comparing it with the prescription ()3.24|) based on 
the analysis of the one cut model. 

First of all, for Nf > 2N Cj i, the first cut is always closed. We will see below that in these 
cases with a closed cut the superpotential can be reproduced by setting Si = by hand in 
the corresponding matrix model, confirming the prescription (|3.24|) for Nf > 2iV Cjl . 

Secondly, the lines for Nf = 3 and U{2) x U(l) correspond to the N Cyl < N f < 2N cA 
part of the prescription (|3.24)1 . There indeed are both an open cut solution and a closed 
cut solution. We will see below that the superpotential of the closed cut solution can be 
reproduced by setting S± = by hand in the corresponding matrix model. On the other 
hand, the superpotential of the open cut solution can be reproduced by not setting Si = 0, 
namely by treating Si a dynamical variable and extremizing W e g with respect to it. In fact 
these two solutions are baryonic branch for a closed cut and non-baryonic branch for an open 
cut. 

Finally, for Nf < N Cj i, the cut is always open, which is also consistent with the prescrip- 
tion (|3.24J) . In this case, the superpotential of the open cut solution can be reproduced by 
extremizing W e & with respect to it, as we will see below. For Nf = N Ci i there should be an 
Si = (cii = a"i) solution for some Zf (corresponding to £7(2) theory with Nf = 2 in the 
r = branch) in the quadratic case, but under the constraint (|4.3jl we cannot obtain that 
solution. 

Below we present resulting exact superpotentials, for all possible breaking patterns. For 
simplicity, we do not take care of phase factor of A which gives rise to the whole number of 
vacua. For details of the calculation, see Appendix IU1 



Definitions: W cl = -§ for 17(1) x U{2) and W cl = -\ for U{2) x 17(1). For simplicity we 
set g — 1 and A = a 2 — ai = —m/g = 1. 

• U(I) x U(2) with N f = 1 



Results 



W low = W d 



2T 



5T 2 



115T 3 245T 4 30501T 5 12349T 6 



+ ••■, t = aI. 



2 



+ 



12 4 64 3 



30 



• U(2) x U(l) with N f = 1 

• U(l) x U(2) with N f = 2 

32T 6 3136T 12 
Wiow = W d + 2T 2 - 6T 4 40T 8 - 192T 10 + • • • , T = A. 

• U(2) x U(l) with N f = 2 

1 fiT 6 224DT 12 
W low = W d - 2T 4 - — 24T 8 - 128T 10 - + • • • , T = A. 

o o 

• U(i) x U(2) with Nf = 3 

19T 2 51T 3 157T 4 5619T 5 33T 6 , s 

W low = W cl + 2T- — + — + — r + ^ r + — + ..., T = A2. 

• U(2) x U(l) with N f = 3: two solutions 

5T 2 39fi c >91T 6 
^low.baryonic = W c l + T 33T 3 - 543T 4 - 10019T 5 + • • • , T = A 3 . 

^low = ^cl + A 3 . 

• U(T) x U(2) with N f = 4 

176T 3 

py low = W cl + 2T- 13T 2 + 138T 4 + 792T 6 

3 

-9288T 8 + 137376T 10 - 2286144T 12 + • • • , T = A. 

• U(2) x U(l) with N f = 4 

40T 3 4928T 6 

W low = W c \ + T — 6T 2 56T 4 - 288T 5 9984T 7 - 63360T 8 

3 3 

1244672T 9 

2782208T 10 - 19009536T 11 • • • , T = A 2 . 

• U(l) x U(2) with N f = 5 

_ 33T 2 1525T 3 3387T 4 314955T 5 

w..w = ^-2r-- - — 

74767T 6 _ i 

■ + •••, T = A5. 



3 

• U(2) x U(l) with N f = 5 



19T 2 154T 3 

W low = W cl + T — + — 132T 4 + 828T 6 + • • • , T = A. 
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4.2 Matrix model computation of superpotential 

In this subsection we compute the superpotential of the system ()4.1|) in the framework 
of lj. If all the Nf flavors have the same mass rrif — —Zf, the effective glueball superpo- 
tential W e ff(Sj) for the pseudo-confining phase with breaking pattern U(N C ) — > YYi=i U(N C ^), 
Eti N c,i = N c is, from (£12), 



W, 



off 



1 A N f ( Nf\ Nf 

to) = -5 E ^ - "f n / + [ N ° ~ T ) W(Ao) + 2 

1=1 ^ ' 



2?rz ( iV c - ^ ) 3 + 2nir S + 2vri ^ 6^ 



n-l 



i=l 



where the periods associated with adjoint and fundamentals are defined by 



Mo M Mo 

n,(^-) = 2 / y(z)cfe, n/^-) = / y{z)dz = - / 



y{z)dz, 



y{z) = VWW^W) ■ 



(4.10) 



For cubic tree level superpotential (j4.1j) . the periods 11^2 (Sj) were computed by explicitly 
evaluating the period integrals by power expansion in [3T], as 



ni _j_ 

2#A 3 #A : 



-2s? + 10sis 2 -5sl) 



+ 2s 2 logA c 



i + log M 

+ 91s?s 2 - 118sis 2 + ^s 2 



+ 



280 4 3484 - 2 2 5272 , 871 , , 

--3-S1 + -3-^2 - 2636s?s* + — Sl sl - —4 ) 



IT, 



:[W(A )-W(a 2 )] + s 2 



2#A 3 5A 3 

+ (2s 2 2 - 10 Sl s 2 + 5s?) 



1 + log 


m 







2si log A 



,'280 , 3484 , 
( -3-sl - — s i4 + 2636s^ a - — 



,, > 5272 , 871 . 

^2 + —4] 



where A = a 2 — ai, Sj = Si/gA 3 , and A = A /A. 

Under the constraint (|4.3jl . the contours defining II x and U/ coincide, so 



;Hi = -n. 



A 



y(z)dz. 



(4.11) 



(4.12) 



/= a i 
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Using this, we can rewrite (j4.10j) as 



WeffOSi, S 2 ) = N CiX W{ ai ) + N Cj2 W{a 2 ) 



N f 



W( ai ) - W(z f ) 



-n x - W(A ) + W(a x ) 



N, 



c.2 



L -n 2 - W(A ) + W(a 2 ) 



- 2ni(N cA + N C)2 )S + 27Tir S + 2mbiSi 



(4.13) 



Here we rearranged the terms taking into account the fact that the periods take the form 
= W(A ) — W(ai) + (quantum correction of order O(Si)), and also the fact that we 



are considering Zf 



ai (thus the second term). The first line corresponds to the 



classical contribution, while the second and third lines correspond to quantum correction. 
Furthermore, we defined iV Cjl = N C) i — Nf/2. 

We would like to extremize this W e s (I4.13j) with respect to S± j2 , and compute the low energy 
superpotential that can be compared with the W\ ow obtained in the previous subsection using 
gauge theory methods. In doing that, one should be careful to the fact that one should treat 
the mass Zf as an external parameter which is independent of Si, 2 although we are imposing 
the constraint ()4.3|) . Zf — — a^(Si,S 2 ). Where is the Zf dependence in ()4.13j) ? Firstly, 
Zf appears explicitly in the second term in (j4.13|) . Therefore, when we differentiate W e g with 
respect to Si t2 , we should exclude this term. Secondly, there is a more implicit dependence 
on Zf in 11/ = — y(z)dz, which we replaced with — TIi/2 using (j4.12j) . If we forget to 
treat Zf as independent of Si, then we get an apparently unwanted, extra contribution as 

Ao dy(z) 



= — ^^-dz — ^--y(z)\ z=Zf . However, this last term actually does not make difference 
because 



dzt 



y(z)\ z=Zf = g yj(z- z f )(z - a+)(z - a 2 )(z - a+) 



0. 



= 2 / 



Therefore what one should do is: i) plug the expression (j4.11|) into ()4.13|1 . ii) solve the equation 
of motion for S\ t2 using (j4.13j) without the second term, and then iii) substitute back the value 
of Si t2 into ()4.13|) . now with the second term included. 

Solving the equation of motion can be done by first writing the Veneziano-Yankielowicz 
term (log and linear terms) as 



N c ,t 



-n x - W(A ) + W{a x ) 



N, 



c.2 



-n 2 - W(A ) + W(a 2 ) 



- 2ni(N cA + N C)2 )S + 2niT S + 27ri6 1 S 1 
gA 3 { iV C)lSl [l - Iog( Sl /A?)] + iV Cj2S2 [l - log( S2 /A^)] + Otf) 
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where 



A 



3JV 0j2 



3iV c ,i 



1 




and then solving the equation of motion perturbatively in Ai 2- In this way, one can straight- 
forwardly reproduce the results obtained in the previous section in the case with the first cut 
open. In the case with the first cut closed, in order to reproduce the results in the previous 
section, one should first set Si = by hand, and then extremize W e s with respect to the 
remaining dynamical variable S 2 . 

Following the procedure above, we checked explicitly that extremizing W eff (S'i, S 2 ) (open 
cut) or W e ff(Si = 0, S 2 ) (closed cut) reproduces the W\ ow up the order presented in the previous 
section, for all breaking patterns for U(3) theory. 

In the above, we concentrated the explicit calculations of effective superpotentials in U (3) 
theory with cubic tree level superpotential. These explicit examples are useful to see that 
the prescription ()3.24j) really works; one can first determine using factorization method when 
we should set Si = by hand, and then explicitly check that the superpotential obtained by 
gauge theory can be reproduced by matrix model. 

However, if one wants only to show the equality of the two effective superpotentials on 
the gauge theory and matrix model sides, one can actually prove it in general cases. In 
Appendix El we prove this equivalence for U(N C ) gauge theory with an degree (k + 1) tree 
level superpotential where k + 1 < iV c . There, we show the following: if there are solutions to 
the factorization problem with some cuts closed, then the superpotential W\ ow of the gauge 
theory can be reproduced by extremizing the glueball superpotential W e s{Si) on the matrix 
model side, after setting the corresponding glueball fields Si to zero by hand. Note that, on 
the matrix side we do not know when we should set Si to zero a priori; we can always set Si 
to zero in matrix model, but that does not necessarily correspond to a physical solution on 
the gauge theory side that solves the factorization constraint. 

5 Conclusion and some remarks 

In this paper, taking M = 1 U (N c ) gauge theory with an adjoint and flavors, we studied the 
on-shell process of passing Nf flavor poles on top of each other on the second sheet through a 
cut onto the first sheet. This corresponds to a continuous transition from the pseudo-confining 
phase with U(N C ) unbroken to the Higgs phase with U (N c — Nf) unbroken (we are focusing 
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on one cut). We confirmed the conjecture of [I] that for Nf < N c the poles can go all the 
way to infinity on the first sheet, while for Nf > N c there is obstruction. There are two types 
of obstructions: the first one is that the cut rotates, catches poles and sends them back to 
the first sheet, while the second one is that the cut closes up before poles reach it. The first 
obstruction occurs for N c < Nf < 2N C whereas the second one occurs for N c < Nf. 

If a cut closes up, the corresponding glueball S vanishes, which means that the U(N C ) 
group is completely broken down. This can happen only by condensation of a charged massless 
degree of freedom, which is missing in the matrix model description of the system. With a 
massless degree of freedom missing in the description, the S = solution should be singular 
in matrix model in some sense. Indeed, we found that the S = solution of the gauge theory 
does not satisfy the equation of motion in matrix model (with an exception of the Nf = N c 
case, where the S = solution does satisfy the equation of motion). How to cure this defect 
of matrix model is simple — the only thing the missing massless degree of freedom does is 
to make S = a solution, so we just set S = by hand in matrix model. We gave a precise 
prescription (J3.24)) when we should do this, i.e., in the baryonic branch for N C j < Nf < 2N c j 
and in the r = N c j non-baryonic branch, and checked it with specific examples. 

The string theory origin of the massless degree of freedom can be conjectured by generaliz- 
ing the argument in |HJ. We argued that it should be the D3-brane wrapping the blown up S 3 , 
along with fundamental strings emanating from it and ending on the noncompact D5-branes 
in the Calabi-Yau geometry. 

Although we checked that the prescription works, the string theory picture of the S = 
solution needs further refinement, which we leave for future research. For example, although 
we argued that some extra degree of freedom makes S = a solution, we do not have 
the precise form of the superpotential including that extra field. It is desirable to derive it 
and show that S = is indeed a solution, as was done in [H] in the case without flavors. 
Furthermore, we saw that there is an on-shell S = solution for Nf = N c . Although this 
solution solves the equation of motion in matrix model, there should be a massless field behind 
the scene. It is interesting to look for the nature of this degree of freedom. It cannot be the 
D3-branes with fundamental strings emanating from it, since for this solution the noncompact 
D5-branes are at finite distance from the collapsed S 3 and the 3-5 strings are massive. Finally, 
we found that the S = solution is in the baryonic branch. It would be interesting to ask if 
one can describe the baryonic branch in the matrix model framework by adding some extra 
degrees of freedom. 

In the following, we study some aspects of the theory, which we could not discuss so 
far. We will discuss generalization to SO(N c ) and USp(2N c ) gauge groups by computing the 
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effective superpotentials with quadratic tree level superpotential. 



5.1 SO(N c ) theory with flavors 

Here we consider the one cut model for SO(N c ) gauge theory with Nf flavors. The tree 
level superpotential of the theory is obtained from M = 2 SQCD by adding the mass for 
the adjoint scalar $ 

W tree = yTr^ + Q^QO/f + Q^/f^' (5.1) 

where / = 1, 2, • • ■ , 2Nf and the symplectic metric Jff and mass matrix for quark rhff> are 
given by 

J = ( ^ Q^j ® l N f xN f , fh = ( ^ g ^ ® diag(mi, • • • , m Nf ) . 

For this simple case the matrix model curve is given by 

y{zf = m\ (z 2 - V) . (5.2) 

This Riemann surface is a double cover of the complex z-plane branched at the roots of y^ 
(that is z = ±2/i). 

The effective superpotential receives contributions from both the sphere and the disk 
amplitudes in the matrix model jB] and the explicit form was given in jT] for U(N C ) gauge 
theory with flavors. Now we apply this procedure to our SO(N c ) gauge theory with flavors 
and it turns out the following expression 



If" \ f 1 2Nf f Ao I 

~2 ( E N a - 2 J j Sr y( z ) dz - 4 E/ y^ dz + 2 {2Nc - 4 - 2Nf) w{Ac 
. 2N f 

+- J2 w ( z i) - ni ( 2N c - 4 - 2N f) s + 2mT 3 + 2m h s i 



eff 



1=1 i=l 

where S = S + 2 Ym=i &i an d z i is the root of 



B(z) = det m(z) = J[ (z 2 - zj) . 



i=i 



Since the curve f!5.2|) is same as the one (j3.1|) of U(N C ) gauge theory, we can use the 
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integral results given there to write down the effective superpotential as 

(N c - 2) 



W, 



eflf 



s 



2 

N f 



E 

7=1,77=0 

s E 

j=l, r/ =l 



+ log 
-log 
-log 



2V m ^A*- 2 " w /detz 



S 1 



2 + 2 



f 1 
2 ~ 2 



m^zf / 25 



2g 



2ff 

m A z\ 



+ 



m A zj 
2S 



-4 1 



2ff 



- f 



We can solve M(z) and T(z) as did for U(N C ) gauge theory. For simplicity we take all 
rj = 0, i.e., all poles at the second sheet. For the J-th block diagonal matrix element of M(z) 
(/ = !,••• , Nj) it is given by 



MAz) 





R(z)—R{qj=—mj) 
z+mj 



R(z)—R(qi=mi) 
z—mj 





where 



R(z) = txia (z — \/ z 2 — . 

Expanding Mz(z) in the series of z we can find (Qf$Q?J ffl + Q f m ff ,Qf) = 2N f S. 
The gauge invariant operator T(z) can be constructed similarly as follows: 



T(z) 



B'(z) 



2B 



-E 



y{qi)zi 



c{z) 2R{z) 
(z) y(z) (z 2 - z 2 ) ' y(z) z y(z) 



+ 



(5.3) 



where 



c(z) 



Tr 



iV f 



z — $ 



For the theory without the quarks, the Konishi anomaly was derived in ^SJElEll- The last 
term in ()5.3|) reflects the action of orientifold. For our example we have 



and 

T(z) 



1 

z 

+- 



Nt 



c{z) = m A (N c - N f ) . 



N f , 

Y^zAzt-^jz*- VJ + 2fl 2 (JV c -2-JV» 

7=1 ^ ' 

TV". 



E ^ 4 - E v z2 - 4 ^ ( z * + 2 ^ 2 ) + 6 ^ - ^/) - 12 » 4 



i=i 



i=i 
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f) 



where for equal mass of flavor, we get (Tr $ 2 ) = Njq (^q — \j q 2 — 4/i 2 j + 2/i 2 (iV c — 2 — N 
and (Tr $ 4 ) = iY^g 4 - N s ^q 2 ~^ 2 q (q 2 + 2/j 2 ) + 6/1 4 (N c -2- N f ). 

Let us assume the mass of flavors are the same and K of them (in this case, 77 = 1) locate 
at the first sheet while the remainder (Nf — K) where 77 = are at the second sheet. Then 
from the effective superpotential, it is ready to extremize this with respect to the glueball 
field S 19 



= log 



A f(JV -2)-JV / ' 



+K log I - v'f + (N f - K) log ( Zf+ V Z f g I . (5.4) 



Or by rescaling the fields S = 2 m A k 2 > % = "a one S e ^ s ^ ne solution and consider for X = 
case 



N f 



This equation is the same as the one in U(N C ) case with N c — > ^p^, so the discussion of 
passing poles will go through without modification and the result is when Nf > Nc ~ 2 , the 
on-shell poles at the second plane cannot pass the cut to reach the first sheet far away from 
the cut. 

By using the condition 1)5.40 . one gets the on-shell effective superpotential 



1 ,„ _ „ , „ 1 2 / „ ,_ L 2S 



Weff.on-ahdi = - (N c - 2 - N f ) S + -m A z/ Nf + (2K - Nf) Wl 



It can be checked that this is the same as \m A (Tr$ 2 ). 

5.2 USp(2N c ) theory with flavors 

For the U Sp(2N c ) gauge theory with Nf flavors we will sketch the discussion because most 

of them is similar to SO(N c ) gauge theory. The tree level superpotential is given by ()5.1)) but 

19 One can easily check that this equation with parameters (N c , Nf, K) is equivalent to the one with param- 
eters (N c — 2r, Nf — 2r, K — r). The equation of motion for glueball field is the same. Since the equation of 
motion for both r-th Higgs branch and (Nf — r)-th Higgs branch is the same, one expects that both branches 

have some relation. By redefinition of £ — > 4m g A = S,Zf — > 2m £ A Zf = z) we get the final relation between 
K Higgs branch and (Nf — K) Higgs branch. 



38 



with J the symplectic metric, and mff the quark mass given by 

J=\ Z 1 <S)In c xn c , fh = ^ ^ q 1 ^ ® diag(mi, • • • ,771^) • 



We parametrize the matrix model curve and the resolvent R{z) as 

y 2 m = W\zf + f{z)=m\ (^ 2 + V), 
R(z) =m A (z- v 7 ^ 2 + V) . 

The effective superpotential is given by 



W eS = S(N c +l) 




2S \ 


m A z] 


m A z 2 j J 


2S 


2S \ 


m A z] 


m A z 2 j j 


2S 




2S A 1 

m A zj / 2 



1 + 



2S 
m A z 2 j 




and T(z) is given by 

T(z) 



b'(z) y((ii) z i 



+ c(z)_ + 2R(z) 



2B(z) ^y( z )(z 2 -z 2 I ) y(z) z y(z)' 



Note the last term (different sign) compared with the SO(N c ) gauge theory. From the solution 
of M(z), we can show that although (Qf($J) ff/ Qf' + Qfm ff/ Q f 'j) = 2N f S ^ 0, we still 
have on-shell relation \m A (Tr<3> 2 ) = W c g. 
The equation of motion is given by 20 



, % - + 4:S\ z f + Jz) + AS 

= log S-^' 1 + Klog[ V \+( N f- K ) log ' ' 



where S = 2m ^ A 2 , % = ■ From this we can read out the following result: when Nf > N c + 1, 
the on-shell poles at the second sheet cannot pass through the cut to reach the first sheet far 
away from the cut. 



20 One can easily check that this equation with parameters (2N C , Nf, K) is equivalent to the one with 
parameters (2N C — 2r, Nf — 2r,K — r). In other words, the equation of motion for glueball field is the same. 
Since the equation of motion for both r-th Higgs branch and (K — r)-th Higgs branch is equivalent to each 
other, one expects that both branches have some relation. 
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Appendix 

A On matrix model curve with Nf(> N c ) flavors 

In this Appendix we prove by strong coupling analysis that matrix model curve corre- 
sponding to U (N c ) supersymmetric gauge theory with Nf flavors is exactly the same as the 
one without flavors when the degree (k + 1) of tree level superpotential W tre e is less than N c 
21 . This was first proved in [32| but the derivation was valid only for the range Nf < N c . 
Then in (THj, the proof was extended to the cases with the range 2N C > Nf > N c . However, 
in [IB], the characteristic function Pn c (x) was defined by Pn c ( x ) = det(x — $), without taking 
into account the possible quantum corrections due to flavors. In consequence, it appeared 
that the matrix model curve is changed by addition of flavors. In this Appendix, we use the 
definition of Pn c ( x ) proposed in Eq. (C.2) of [1: 

P H .(x) = Axp (- ± fj + A 2 "-""|gexp (± fj , (A.1) 

which incorporates quantum corrections and reduces to PnX x ) = det(x — $) for Nf = 0, and 
see clearly that the matrix model curve is not changed, even when the number of flavors is 
more than N c . Since Pn c {x) is a polynomial in x, (lA.ljl can be used to express U r with r > N c 
in terms of U r with r < N c by imposing the vanishing of the negative power terms in x. 

Assuming that the unbroken gauge group at low energy is U(l) n with n < k, the factor- 
ization form of Seiberg-Witten curve can be written as, 

P 2 Nc (x)-4A 2 N <- N tB(x) = H% c _ n {x)F 2n {x). 

21 The generalized Konishi anomaly equation of R(z) given in l|2.6|l is same with or without flavors, so the 
form of the solution is the same for gauge theory with or without flavor. In this Appendix we use another 
method to prove this result. 
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The effective superpotential with this double root constraint can be written as follows 22 . 

k 

W eS = g r U r+ i 



r=Q 
N c -n 



i=l 



Pv c (x)-2e,A 



N c 



B(x) 



-dx + Bi 



P Nr (x)-2e i K 



n ^a/B(x) 



x-pi 

The equations of motion for Bi and Pi are given as follows respectively: 
P N (x) -2e,A 



S -Pi) 



-dx 







B(x) 



dx, = 2R- 



Pat. x -2e,-A 



N c -4- 



B{x) 



dx. 



(x-p^ 2 J (x-pi) 3 

Assuming that the factorization form does not have any triple or higher roots, we obtain 
Bi = at the level of equation of motion. Next we consider the equation of motion for U r : 

N c -n 



= 9r-l + 



i=l J 



x r x r \ a; 4 



i=i 



x-_p. 



dP N 



where we used Bi = and (jA.lj) to evaluate du 
and sum over r. 



Now, as in jS], we multiply this by z r 1 



A',, 



x — Z 14 — ' X — Pi 

1=1 ^ l 



2x 



x fc / ^ 



N c -n 



X — Z 



-exp 



-dx. 



fe=i 



7 X-Pi 



Defining the polynomial Q(x) in terms of 

Nc-n 



E 



g(x) 



x -p» H Nc - n {x) 
and also using (jA.ljl and factorization form, we obtain 



(A.2) 



W\z) 



P 



Q(x) 



x - z H Nc _ n (x 
Q(x)^/F 2n (x) 



-dx + 



P 



N c 



Q(x) 



x - zH Nc Jx 



-dx + 



dx 



x — z 



dx. 



x — z 



This is nothing but (2.37) in j^j. Since W'(z) is a polynomial of degree /c, the Q should be 
a polynomial of degree (k — n). Therefore, we conclude that the matrix model curve is not 
changed by addition of flavors: 



y 2 m = F 2n (x)Qi_ n (x) = W' k {xY + 0(x 



Jfc-i^ 



(A.3) 



22 If we want to generalize this proof to more general cases in which k + 1 is greater than and equals to N c , 
we have to take care more constraints like Appendix A in 8;, which should be straightforward. 
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B Equivalence between W\ ow and W e s((Si)) with flavors 



In this Appendix we prove the equivalence W\ ow in U(N C ) gauge theory with W e s((Si)) in 
corresponding dual geometry when some of the branch cuts on the Riemann surface are closed 
and the degree (k + 1) of the tree level superpotential Wtree is less than N c . This was first 
proved in [32], however, the proof was only applicable in the Nf < N c cases. Especially, the 
field theory analysis in [32] did not work for iV c < Nf < 2N C cases. Furthermore, as we saw 
in the main text, for some particular choices of Zj (position of the flavor poles), extra double 
roots appear in the factorization problem. In section HJ we dealt with £7(3) with cubic tree 
level superpotential and saw the equivalence of two effective superpotentials for such special 
situations. To include these cases we are interested in the Riemann surface that has some 
closed branch cuts. Therefore our proof is applicable for U (N c ) gauge theories with W tre e of 
degree k + 1 (< N c ) in which some of branch cuts are closed and number of flavors is in the 
range iV c < Nf < 2N C . In addition, we restrict our discussion to the Coulomb phase. 

In the discussion below, we follow the strategy developed by Cachazo and Vafa in j3J)] and 
use (IA.1J) as the definition of Pn c (x). We have only to show the two relations: 

W low (g r , Zl ,A)\ A ^ = W^((Si))\ A ^ (B.l) 
dW low (g r , Zl ,A) dWesdSj)) 

dX = — ' (R2) 

the equivalence of two effective superpotentials in the classical limit and that of the derivatives 
of the superpotentials with respect to A. 

B.l Field theory analysis 

Let k be the order of W( Tee and n(< k) be the number of U(l) at low energy. Since we 
are interested in cases with degenerate branch cuts, let us consider the following factorization 
form 23 : 

P Nc (x) 2 - 4A 2N <- N tB(x) = F 2n (x) [Q k - n (x)H Nc _ k (x)~\ * = F 2n (x) [H Nc „ n (x)} 2 , (B.3) 
W'(x) 2 + / fc _i(x) = F 2n (x)Q k „ n (x) 2 . 

If k equals to n, all the branch cuts in F 2k {x) are open. The low energy effective superpotential 
23 In the computation below, we will use relation 1|A.2|) and put g k +i = 1- 
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is given by 



n+l 



r=l 



I 



1=1 



U [ PnM - 2e i A N ^y/B(p i ) ) + Q~ f P Nc (Pi) - 2e i h N ^ \j B{ Pi ) 







where I = N — n and Pn c (x) is defined by (C.3) or (C.4) in |T], 
P Ne (x) = (det(x-$)) + 



A 2N c - Nf B{x)_ exp 



X 



x i 
i=l 



(B.4) 



The second term is specific to the Nf > N c case, representing quantum correction. Define 

K(x) 



A 2N c ~Nf B\ x ) 

A c 1 — — exp 



x 



N c 



The first term in (jB.4|) can be represented as (det(x — $)) = X]fc=o % Nc ~ k Sk- The relation 
between C/j's and s^'s are given by the ordinary Newton relation, ks^ + Yl r =i r ^r s k-s — 0. 
From the variations of W\ ow with respect to Pi and Q iy we conclude that Q, L = at the level 
of the equation of motion. In addition, the variation of W\ OVf with respect to U r leads to 



i=i 



j=l j=0 

/ 



j=l 



i=l j=0 



i=l 



B(P: 
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2N r -N 
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By using these relations, let us compute W' cl 

r=l 

Nr I Nr . I N c I 



r=— oo i=l j=0 i=l r=l i=l 

E ^TT^ 1 " " 5 E " •) " E E i.ft M*" 1 

i=l «=1 r=l i=l 

£ - E - - E + - E 

i=l ^ i=l ^* j=l i=l 



r=l i=l 

-2> 



where we dropped 0(x ). The fifth term above can be written as 

N c I N c I I 

-EE L <^) xr ~ 1 = - E E L ^^) ;r " 1 + E- L ^) + ( 

r=l i=l r=—oo i=l i=l 



X 2 ) 



After some manipulation with the factorization form (jB.3|) . we obtain a relation 

Qk-n{x)VF 2 n{x) = -~ ~ -— exp > 

H Nc - k H Nc _ k x N 4 



X 1 
1=1 



_ Pnc{x) _ 2K{x) 

HNc-k HNc-k 

Substituting this relation into (jB.5|) we obtain 

i 



0(af 



W' A = Q k -n(x) v /F 2n (x)- Y] - [LiP N M ~ 2LiK(pi)] + V ^-L t 

— ' ~^ x — Pi 



i=l 
N c I 

~ E Y^LiGripi^ + Oix-*) 
r=—oo i=l 

Finally let us use the following relation 24 ; 



l , s n c i 

y^K[X)_ h _ y, y- Li g M)x r-l = 0{x -2y ( B .6) 

^ X-Pi , 
i=l r=—oo i=l 



24 For simplicity, we ignore Ei=i ^i- P rove the relation, let us consider the circle integral over C . Until 
now, since we assumed pi < x, the point pi should be included in the contour C . Multiplying pr fc > and 
taking the circle integral, we can pick up the coefficient Ck-i of x k ~ 1 . In addition, if we denote a polynomial 
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After all, by squaring W' ch we have 

i 

Qk-n(x) 2 F 2n (x) = W%(x) + 2 ~ [LiPNXVi) ' 2L t K( Pl )] x^ 1 + 0(x k ~ 2 5 

i=i 

1 1 



»=i 



On the other hand the variation of with respect to A is given by 

Wow T „/ x 1 T „ , x &fc 



dlogA 2 ^"^ ^ 2 

This is one of the main results for our proof. In the dual geometry analysis below, we will see 
the similar relation. 

In the classical limit, we have only to consider the expectation value of In our assump- 
tion, gauge symmetry breaks as U(N C ) — > U(N Cji ) we have Tr $ = Y17=i Nc^i- Therefore 
in the classical limit W\ ow behaves as 

n 

W cl = J2 N c,iW( ai ). (B.8) 

i 

By comparing (|B.7|) and (jB.8|) and the similar result which we will see in the dual geometry 
analysis below, we will show ()B.1|) and (jB.2|) . Let us move to the dual geometry analysis. 



B.2 Dual geometry analysis with some closed branch cuts 

As we have already seen in the main text, a solution with (Si) = appears for some special 
choice of zi, the position of flavor poles. In our present proof, however, we put some of Si 
to be zero from the beginning, without specifying zj. More precisely, what we prove in this 



M = A 2Nc Nf -pi^-exp (X^i §0 = Sfc-oo a 3 x ^ we can obtain following relation, 



[M] + 



§1 + k fa jX -^ <=► ^=&(,) + Ea^) 



J + 3=0 3=0 

where right hand side means circle integral of left hand side. Thus, we obtain the Ck-i as 



c fc _i=V/ ^ dx + V/ ^ dx=~y^(f ajx 3 k V -^rdx + V a 3 p\ k = 

where we used that around x = 0, x ^ p = — YlnLi ^ rr - Therefore the left hand side of (IB. 61) can be written 



as E^L-oo c 3 xi = £,•=?-«, c J xj = °( x 2 )- 



45 



Appendix is as follows: for a given choice of zj, if there exists a solution to the factorization 
problem with some of (Si) vanishing, then we can construct a dual geometry which gives the 
same low energy effective superpotential as the one given by the solution to the factorization 
problem, by setting some of Si to zero from the beginning. Therefore, this analysis does 
not tell us when a solution with (Si) = appears. To know that within the matrix model 
formalism, we have to go back to string theory and consider an explanation such as the one 
given in 

Now let us start our proof. Again, let k be the degree of tree level superpotential W{ ree (x) 
and n be the number of U(l) at low energy. To realize this situation, we need to consider 
that (k — n) branch cuts on the Riemann surface should be closed, which corresponds to 
(Si) = 0. Here, there is one important thing: As we know from the expansion of W e Q in terms 
of A (e.g. see (j4.13|) ). we cannot obtain any solutions with (Si) = if we assume that Si is 
dynamical and solve its equation of motion. Therefore to realize the situation with vanishing 
(Si), we must put Si = at the off-shell level by hand. With this in mind, let us study dual 
geometry which corresponds to the gauge theories above. In the field theory, we assumed that 
the Riemann surface had (k — n) closed branch cuts. Thus, in this dual geometry analysis 
we must assume that at off-shell level, (k — n) S^s must be zero. For convenience, we assume 
that first n S^s are non-zero and the remaining (k — n) vanish, 

Si 7^ 0, i = l,---n, Si = 0, i = n + l,---k. 

Therefore the Riemann surface can be written as 

V 2 = F 2n {x)Q 2 k _ n = W'(xf + + ■■■ 

The effective superpotential in dual geometry corresponding to U(N C ) gauge theory with Nf 
flavors was given in [I] (See also (j2.12j) ) and in the classical limit it behaves as 25 

n 

W c z\ cl = Y, N c,iW(ai). (B.9) 

i=l 

As discussed in the previous section, existence of flavors does not change the Riemann 
surface y(x). In other words, Riemann surface is not singular at xj (roots of B(x)), 

& y(x)dx = 0, y(x) = ^W'(x) 2 + b k ^x k - 1 H . 

J Xj 

25 Remember that in this Appendix we are assuming only Coulomb branch. For the Higgs branch, see (7.11) 
and (7.12) in Q]. 
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Therefore as in [35] . by deforming contours of all S^s and evaluating the residue at infinity 
on the first sheet, we obtain the following relation, 

n k j 

^2 Si = ^2 Si = -jh-u 

i=l i=l 

where we used Y^H= n +i &i — 0- 

With this relation in mind, next we consider the variation of W e s with respect to Si'. 

— = 0, i = l---n. (B.10) 

Solving these equations, we obtain the expectation values, (Si). Of course, these vacuum ex- 
pectation values depend on A, g r and zj. Thus when we evaluate the variation of W e ff ((»%), A) 
with respect to A, we have to pay attention to implicit dependence on A. However the implicit 
dependence does not contribute because of the equation of motion ()B.10|) : 

dW eS ((Sj),A) = ^d(Sj) dW^jS^A) dW^jS^A) 

dA 2-? OA d(Si) OA 

i=i 1 ' 

dW eS ((Sj),A) 
OA 

On the other hand, explicit dependence on A can be easily obtained by monodromy analysis. 
Here let us recall the fact that the presence of fundamentals does not change the Riemann 
surface. In fact, looking at (|2.12jl we can read off the dependence from the term 2ttitq = 



log 



( b l a 2 



V A 



->.W 



d l ogA 2N c -N f — 

To finish our proof, we have to pay attention to fk-i{x), on-shell. Namely putting (Si) 
into fk-i(x) what kind of property does it have? To see it, let us consider change of variables 
from SiS to b^s. As discussed in jH] the Jacobian of the change is non-singular if < j < k — 2, 

OS, 1 f , x j 
ax- 



dbj 8m J Ai ^ W '(x) 2 + f(x)' 

In our present case, since only n of k S^s are dynamical variable, we use 6j, i — 0, • • • n — 1 in 
a function f k _i(x) = ^biX 1 as new variables, instead of Si's. As discussed in jSHlE21IIl, by 
using Abel's theorem, the equation of motions for fej's is interpreted as an existence condition 
of a meromorphic function that has an A^-th order pole at infinity on the first sheet and an 
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(iV c — Nf)-th order zero at infinity on the second sheet of £ and a first order zero at qi. For 
a theory with Nf < 2N C , such a function can be constructed as follows [331 IT]: 

tP( x ) = p Nc ( x ) + yJP^Jx) -4A 2N c~ N fB{x). 

For this function to be single valued on the matrix model curve y(x), the following condition 
must be satisfied, 

P 2 » -4A 2N <- N fB(x) = F 2n (x)H 2 Nc _ n (x) 
W\x) 2 + f(x) = F 2n {x)Q 2 k _ n (x) 

This is exactly the same as the factorization form we already see in the field theory analysis. 
Therefore the value b^-i of on-shell matrix model curve in dual theory is the same one for 
field theory analysis. Comparing two results, (|B.7j) and (|B.8j) with corresponding results for 
the dual geometry analysis, ()B.9|) and (jB.ll|) we have shown the equivalence between these 
two descriptions of effective superpotentials. 

C Computation of superpotential — gauge theory side 

In this Appendix, we demonstrate the factorization method used in subsection 14.11 to 
compute the low energy superpotential, taking the Nf = 4 case as an example. Therefore 
there are two kinds of solutions for the factorization problem ()4.5j) and ()4.6|) . 

• The breaking pattern U(2) x U(l) 

The first kind of solution for the factorization problem is given by 

A = 0, B = 2a, C = a 2 - 4A 2 , D = 0, c = 0, b = 0. 

In the classical limit A — > 0, we can see the characteristic function goes as P^{x) — > x 2 (x + a), 
which means that the breaking pattern is U(2) x U(l). Note that since we are assuming 
rrif = 0, the notation " ~ " should be used for the gauge group that corresponds to the cut 
near the critical point at x — 0. Inserting these solutions into ()4.7)1 we obtain one constraint, 

A 2 = a 2 + 8A 2 . 

We can easily represent a as a Taylor expansion of A: 

a = -1+4T + 8T 2 + 32T 3 + 160T 4 + 896T 5 + 5376T 6 + --- , 
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where we put A = 1 and defined T = A 2 . 26 The coefficients of P${x) are related to the 
Casimirs Uj = }(Tr[$ J ']) as follows. For N c = 3, JV> = 4, (1A~TJ) reads 



P a {x) = x 3 exp (" £ §) + A 5 ^ exp (^ ^ 
\ 7=1 / \ 7=1 ■ / 



3 

U l , A2 \ „ , / r 7" I T T T T U l A2; 



= x d - C/i^ + ( -U 2 + ^- + A z jx + \^-U 3 + UtU 2 - -± - A'Uxj +■■■ . 
Comparing the coefficients, we obtain 

2 3 

Ui = -a, U 2 = -b+ °- + A 2 , U 3 = -c + ab-^--aA 2 . 

Furthermore, one can compute ai j2 from ()4.7|) . Plugging all these into ()4.9|) . we finally obtain 

4DT 3 4Q28T 6 1 

^ low = W cl + T-6T 2 -^-56T 4 - 288T 5 -^^ + ---, T = A 2 , W cl = --. 

6 6 O 



The breaking pattern U(l) x £7(2) 



The other kind of solution for the factorization problem is given by 

A= -(-a-27?A), B = a-2r]A, C = -(a + 2r/A) 2 , 
£1 = 0, c = 0, 6=^(a + 2r/A) 2 

where 77 = ±1. These solutions correspond to the breaking pattern U(l)xU (2) in the classical 
limit. Inserting these solutions into ()4.7|) we obtain one constraint, 

A 2 = ^(a 2 -20ar/A + 4A 2 ). 

Again, let us represent a as a Taylor series of A: 

a = -2 + 10T-24T 2 + 144T 4 - 1728T 6 + -- - , 

where we put A = 1, rj = 1 and defined T = A. Doing the same way as previous breaking 
pattern, we can compute the effective superpotential as 

176T 3 1 

W low = W cl + 2T - 13T 2 + — 138T 4 + 792T 6 + ■ • • , T = A, W d = --. 

o o 

The other cases with iVj = 1, 2, 3 and 5 can be done analogously. 

26 If we take care of a phase factor of A, we will obtain the effective superpotentials corresponding to each 
vacuum. However in our present calculation, we want to check whether the effective superpotentials of two 
method, field theory and dual geometry, agree with each other. Therefore, we have only to pay attention to 
the coefficients in Wi ow , neglecting the phase factor. 
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